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Abstract 

We calculate the two-photon exchange (TPE) corrections in the region where the kinematical 
variables describing the elastic ep scattering are moderately large s ~ — t ~ —it 3> A 2 . For this 
kinematics we use the QCD factorization approach formulated in the framework of the soft-collinear 
effective theory (SCET). Such technique allows us to develop a description for the soft spectator 
scattering contribution which is found to be important in the region of moderately large scales. 

Together with the hard spectator contribution we present the complete factorization formulas 
for the TPE amplitudes at the leading power and leading logarithmic accuracy. The momentum 
region where both photons are hard is described by only one new nonperturbative SCET form factor. 
It turns out that the same form factor also arises for wide-angle Compton scattering which is also 
described in the framework of the SCET approach. This allows us to estimate the soft spectator 
contribution associated with the hard photons in a model independent way. 

The main unknown in our description of the TPE contribution is related with the configuration 
where one photon is soft. The nonperturbative dynamics in this case is described by two unknown 
SCET amplitudes. We use a simple model in order to estimate their contribution. 

The formalism is then applied to a phenomenological analysis of existing data for the reduced 
cross section as well as for the transverse and longitudinal polarization observables. 



Introduction 

The electromagnetic form factors of the nucleon, have been explored extensively during the past 50 years 
with ever increasing accuracy. The tool to extract the electromagnetic form factors is provided by the 
one-photon (I7) exchange approximation to elastic electron-nucleon scattering. Precision measurements 
of the proton electric to magnetic form factor ratio at larger Q 2 using polarization experiments [l]-[6] 
have revealed significant discrepancies in recent years with unpolarized experiments using the Rosenbluth 
technique I7H9], when analyzing both within the 17-exchange framework. This discrepancy between po- 
larized and unpolarized measurements has generated a lot of activity, both theoretical and experimental, 
over the past decade. The most plausible solution of this problem is the correct calculation of the ra- 



diative corrections (RC) to elastic lepton-proton scattering. The well-known calculations in Ref. 10 
were performed in the so-called soft photon approximation which allows one to compute the two-photon 
exchange (TPE) diagrams. This calculation has been recently reconsidered and improved in a series 
of works |llf|17| within a hadronic framework. A more detailed review of the subject can be fond in 
Refs. 18 -20]. However with increasing energy, calculations using hadronic degrees of freedom become 
less and less reliable. At large energy s and momentum transfer t = —Q 2 one has to properly take into 



1 On leave of absence from St. Petersburg Nuclear Physics Institute, 188350, Gatchina, Russia 



account the interactions of multiple hard photons with constituents inside the proton. Such dynamics 



was not considered in the pioneering papers 10 which have been published before the formulation of the 



underlying microscopic theory of the strong interactions QCD. 

A consistent description of the hard and soft QCD dynamics can be carried out using the factorization 
approach which was applied for the analysis of many exclusive reactions, see e.g. 21-23 . However even 
this consideration may be sometimes complicated and challenging. The well known example is the 
description of the nucleon form factors (FFs) at large momentum transfer. For many years a proper 
theoretical framework for this regime has been the subject of theoretical debates. There are two different 
points of view on the description of QCD dynamics at large Q 2 . They are related with the different 
mechanisms of the underlying scattering which we describe as a hard and soft spectator scattering. The 
hard spectator scattering mechanism was studied long time ago, see e.g. 21-23 and references therein, 
and results in the well known factorization formula for the FF F\ . On the other hand the soft spectator 
scattering has only been estimated using various phenomcnological approaches such as light-cone wave 
functions and QCD sum rules [24| - p30] . Such estimates show that at moderate values of Q 2 < 10 GeV 2 a 
dominant contribution to the nucleon FFs originates from the soft overlap mechanism. 

Recently we suggested the factorization formula for the nucleon FFs [3l] which systematically includes 
both contributions. We used the soft collinear effective theory (SCET) framework in order to describe the 
factorization of the soft spectator scattering contribution. In this approach the factorization of the soft 
spectator scattering contribution at moderate values of Q 2 can be defined in terms of SCET FFs which 
can be rigorously defined in the intermediate effective theory SCET-I. If the hard-collinear virtualitics 
are relatively small then the further factorization can not be performed and the SCET FFs must be 
considered as nonperturbative functions. 

In the present work we follow the same approach in order to describe the factorization of the short 
and long distances for the two-photon exchange (TPE) contribution shown in Ihs in Fig[l] This allows us 
to perform an unambiguous and consistent separation of the different regions in the QED loops in Fig{T] 
associated with hard and soft configurations of the photons. We restrict our considerations to the region 
where the Mandelstam variables are much larger then the typical QCD scale s ~ — t ~ — u A 2 . In 





Figure 1: The factorization of the TPE contribution in elastic ep scattering. The crossed box diagram 
on the Ihs is not shown for simplicity. The possible attachments of the soft photon 7 S to the lepton lines 
on the rhs are shown by crosses. 



full analogy with the nucleon FFs, the leading power behavior of the TPE amplitudes can be described 
by two different configurations associated with the soft and hard spectator contributions. Schematically, 
the structure of the leading contributions is shown in Figjl] The first two graphs on the rhs describe the 
soft spectator contribution and the third one corresponds to the hard spectator configuration. The hard 
subprocesses, which include one or two hard photons arc shown as blobs with the symbol H. The hard 
spectator configuration has already been studied in Ref. |32||33] . The analysis of the soft spectator terms 
is new. The soft QCD dynamics for the corresponding diagram is described by a SCET FF denoted by T\ 
and amplitudes <?i 3 which are defined in SCET-I. At leading order in the QCD coupling, we obtain these 
new SCET amplitudes which do not appear in the factorization formulas for the nucleon FFs. It turns 
out that J-\ can be fixed from the wide-angle Compton scattering using the universality of its definition 
in SCET-I. The amplitudes 171,3 involve matrix element which cannot be related to known objects, and at 
present can only be estimated within some model approach. As a first step we consider an estimate using 
the effective theory with hadronic degrees of freedom. In this way the TPE contribution is completely 
defined. We then perform phenomcnological studies of the TPE effects and compare our results with 
existing experimental data. 

Our paper is organized as follows. Section[l]is devoted to the general properties of elastic ep scattering. 
We specify notations and kinematics, discuss the general properties of the amplitudes, as well as the 
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structure of the reduced cross section and asymmetries. In section [2] we derive the leading order SCET-I 
factorization formula for the soft spectator scattering contribution in the TPE amplitudes. In the next 
section we perform the matching and compute the one-loop, leading order in a s hard coefficient functions. 
The section [4] is devoted to the extraction of the second unknown FF T\ from the data for wide-angle 
Compton scattering. In section [5] we discuss the the SCET amplitudes 51,3 which describe the TPE 
contribution when one of the photons is soft. We use a simple hadronic model in order to estimate this 
contribution. 

In section [6] we use the obtained result for the phenomenological analysis and estimate the effect of 
the TPE contribution for different observables. The summary of our obtained results is presented in 
section [7] In Appendices A-D we present more details of some calculations. 



1 Elastic lepton-nucleon scattering at large Q 2 

Let us briefly review the main definitions and some results for the process e(k) +p(p) — > e{k') +p(p'). In 
order to describe this process we introduce following notations 



P = l(p + p>), K= l -{k + k'), q=p'-p = k-k', 



and define the Mandelstam variables 



s = (p + kf, t = q z = -Q z , u={p- k'f, t 



Am 2 



where m is the nucleon mass. For further use, we introduce two more convenient variables 



1 + 2(1 + t) tan 2 



(s - u) 2 + t(4m 2 - t) 



v = (K ■ P) 



(s - u) 2 - t(4m 2 - t) ' 
s — u 



< e < 1, 



(1) 



(2) 



(3) 



(4) 



where 9 is the electron Lab scattering angle. One can choose any two independent variables for the 
description of the physical amplitudes of the process. It is customary to use the variables Q 2 and e for a 
description of the cross sections and related observables. 



Then the general parametrization of the ep— scattering amplitude reads 34 

(p' , k' out\ in k,p) = i(2tt) S(p + k — p — k) A ep , 

where 



A ep = —u{k')ru{k) N(p r ) 



rG M (e, Q 2 ) F 2 (e, Q 2 ) + — $ F 3 (e, Q 2 ) 



N{p), 



The canonical one-photon exchange part is given by the well known formula 

„2 



AZ p = —u(k')^u(k) Nip') 



1»G M {Q 2 ) F 2 {Q 2 ) 

m 



N(p), 



(5) 



(0) 



(7) 



The following difference can be considered as definition of the TPE correction arising from the QED 
next-to-leading corrections 



A - Al p ^ Ai; = — 2 u{k')ru{k) 



N{p') 



^SG M (e,Q 2 ) SF 2 ie,Q 2 ) 

m 



$ Fsie, Q 2 ) 



Nip), 



with 



5G M ie, Q 2 ) = G M {e, Q 2 ) - G M (Q 2 ), SF 2 ie, Q 2 ) = F 2 (s, Q 2 ) - F 2 (Q 2 



(8) 
(9) 

(10) 
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The amplitudes SGm, 8F2 and F 3 obtain different contributions from all diagrams associated with the 
QED radiative corrections to elastic ep-scattering. In the present paper we only consider the calculation 
of the leading in 1/Q 2 corrections arising from the TPE contribution, see Figjl] 

Some useful constraints on the behavior of these amplitudes can be established from the consideration 
of their analytical properties. For that purpose, following 16 we introduce the functions: 



G 1 (v,Q 2 ) = 6G E (e,Q 2 ) 



r^ 3 (e,Q 2 ), 



(11) 



G 2 (i/, Q 2 ) = 5Gm(s, Q 2 ) + -^F 3 (e, Q 2 ), 

771 

G 3 (v,Q 2 ) = F 3 (e,Q 2 ), 



(12) 
(13) 

where we use the shorthand notation e = e(v,Q 2 ) and defined SGe = SGm — (1 + t)SF2. 

An analysis of i-channel helicity amplitudes for the ep — ¥ ep process shows that in the Regge limit 
,s ^ 00, Q 2 / s ^ which is equivalent to ^ — >■ 00, Q 2 jv — > 0, the functions Gi vanish 

lim Gi{u,Q 2 ) = 0. (14) 

v—¥oo 

This higher energy behavior allows one to write down unsubtracted dispersion relations for the amplitudes 

,ImG;(i/,Q 2 ) 



d as 16 



dv' 



.~ v' —v 



(15) 



Regge limit can easily be translated into boundary condition for the practically important variable e, and 
corresponds to the limit e 1 at fixed Q 2 . From the Eqs.(ll -13) one obtains: 



lim 



6G EM (s,Q 2 ) + ^F 3 (e,Q 2 ) 



= 0, 



(16) 



lim F 3 (e,Q 2 ) = 0. 



(17) 



Eqs.( 16|17 ) imply that for certain observables the TPE corrections must vanish in the forward limit. 

The analytical expressions for various observables are well known in the literature, see e.g. [34 1 - For 
the convenience of the reader we provide some of them here. The unpolarized cross section reads 



da 

C^Lab 



dam t 
rff7 La b e(l + t) 



(18) 



with the structureless part 



daws 
dh 



4a 2 



E '3 



Lab 



(19) 



The variables E and E' in Eq.(19) denote the incoming and outgoing electron Lab energies. The elastic 
contribution to the reduced cross an section is given by: 



ai{e,Q) = G 



: G 2 E + 2G M Re 



SG 



M 



e^F 3 



2-G E Re 

T 



SG f 



m z 



(20) 



The elastic contribution to the polarization observables measured in the recoil polarization experiments 
reads 



2e(l-e) 



^G e Gm + G E Re 8Gm + Gm Rc 



-A 



)} 



(21) 



<yi.Pi 



e 2 <G 2 M 



2G M Re ( SG 



1 + e m 2 



F? 



(22) 



where P;.t correspond to the recoil proton polarization along or perpendicular to its momentum, re- 
spectively. Besides the elastic contributions, shown in Eqs.( 20p2 |, the observables also include the 
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contribution from the inelastic processes corresponding with the emission of soft photons which provide 
the cancellation of IR-divergent terms in the elastic next-to-leading QED amplitudes. 

Our task is to compute the TPE amplitudes in the limit of large Q 2 ^> A 2 and, if possible, for 
arbitrary values of e. However such task includes the analysis of different kinematical regions which can 
be associated with the different underlying QCD dynamics. Therefore we will split the e-interval into 
three regions which can be described as follows. 

The forward limit e *~ 1 which can be associated with the Regge limit, in which s — > oo and Q 2 is 



fixed. In this case as one can conclude from Eqs.( 16|17 l and Eq.(20) the TPE corrections to the reduced 
cross section are vanishing. 

The second situation is associated with the backward scattering where s ~ Q 2 ^S> \u\ ~ A 2 and 
corresponds to the region of small e « 1. From the unpolarized data we expect that in this case the 
effect of TPE corrections are largest. However, the development of a theoretical approach in systematic 
way for this region is a difficult task. 

The last region is described by kinematics where all Mandelstam variables are large and of the same 
order: 



f > A 2 . 



(23) 



This region can be associated also with the wide-angle scattering. In this situation one can try to compute 
the TPE amplitudes by performing an expansion with respect to the large scale Q 2 with a fixed ratio 
Q 2 1 s. Below we are going to realize this idea using the QCD factorization approach. In this case the 



values of e are restricted to some interval e n 



< e < e n 



where we define the boundaries e min and 



c max from the phenomenological criteria that the minimal value of |u| is given by |u m ; n | =2.5 GeV 2 , 
which guarantees some suppression of subleading power corrections, see Table [I] The upper boundary 



Q 2 ,GeV 2 


3 


4 


5 


6 


£min 


0.60 


0.52 


0.45 


0.42 



can be formally defined by the relevance of the Regge dynamics. For simplicity we will not introduce 
e max assuming an extrapolation to the value e — 1. 

The idea to apply the QCD factorization approach for the wide-angle region has already been used 
in 32 33 where the hard spectator scattering contribution was computed at leading order. Because both 



photons in this case are hard one needs only one gluon exchange as shown in Fig (2] The nonperturbative 
input is described by the nucleon distribution amplitudes (DAs), see the details in Refs. 32 33 . This 





+ 



Figure 2: Reduced diagram describing the hard spectator scattering contribution and the leading order 
diagram with one gluon exchange. 



calculation allows us to estimate the asymptotic behavior of TPE amplitudes at large Q 2 

A^ 

O 6 ' 



5G 



M 



;SF 3 



A 4 
Q 4 ' 



5F 2 



(24) 



The behavior of the amplitudes SGm and 6F3 is similar to the Dirac FF F\ ~ A 4 /Q 4 and can also be 
described within the collinear factorization approach. The helicity flip amplitudes 8F2 is suppressed by a 
power Q~ 2 similar to the corresponding Pauli FF F2 and can not be described by collinear factorization 
due to the end-point divergencies. Therefore qualitatively, upon neglecting the logarithmic structure, the 
situation is quite similar to the nucleon FFs. 
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On the other hand the analysis of the soft spectator scattering for the nucleon FFs shows that 



these terms are not suppressed by inverse powers of Q 35 36 and therefore can also provide sizable 
contributions especially in the region of intermediate Q where QA ~ m 2 . In Refs. [3lj[37] we investigated 
the soft spectator contribution and suggested the generalization of the factorization which includes both 
hard and soft spectator terms. Taking into account that the TPE dynamics is quite similar to the FF case 
one may expect that the same situation is relevant for this case too. Then the soft spectator scattering 
must also be included into the consideration when calculating the TPE amplitudes. 

Such an attempt has been done in 38 where the diagrams in Fig{2]have been evaluted in the framework 
of the so-called handbag approach [54] - In this model the overlap of the hadronic states are described 




Figure 3: Diagrams describing the soft spectator scattering. Both photons interact with the same quark. 

by the generalized parton distribution (GPD) introduced as matrix element of the light-cone twist-2 
matrix element. The GPD which arises here at large Q is considered as a natural generalization of 
the matrix element which appears in the description of the deeply virtual Compton scattering (DVCS) 
at small momentum transfer Q ~ A. However in DVCS kinematics, GPDs describe the soft dynamics 
of the system constructed from collinear partons with a small invariant mass (small Q ~ A ). But 
at large momentum transfer ( (J > A ) one is faced with the different situation which may be better 
associated with the dynamics of jets. In that case the factorization of the soft dynamics is described 

In Ref. 



within the framwork of soft collinear effective theory (SCET) [39 -44 



31 we carried out a 



detailed consideration of the nucleon FFs within the SCET formalism. In the present paper we will show 
that the situation with TPE is quite similar, allowing for a QCD factorization in a similar way. 



2 Soft spectator interaction for TPE in SCET 

The main feature of the SCET factorization is the presence of two large scales: hard ~ Q 2 and hard- 
collinear ~ AQ. At fist step one integrates out the hard modes and passes from the QCD to the effective 
theory SCET-I. If the value of the hard-collinear scale is not large, for instance, AQ ~ m 2 then one 
can not follow further with the factorization. We shall define such a region of Q 2 ~ 2.5 — 16 GeV 2 as 
intermediate. At present, this region covers all existing and planned experiments of nucleon FFs. 

Our task is to construct SCET-I factorization for the TPE amplitudes in this region. Such factorization 
scheme includes the description of the soft spectator contribution within SCET-I and also, if it is necessary, 
the hard spectator contributions. This terms can be considered as part of the SCET-II factorization but 
because their factorization scheme is different they, for instance, do not have the Sudakov suppression, 
they can be considered separately. 

Technically the SCET factorization is constructed by matching of QCD diagrams onto appropriate 
operators in SCET-I effective theory. The matrix elements of the SCET-I operators describe the low scale 
processes which can not be computed in perturbative QCD. The simple power counting arguments allow 
one to expect that the leading order SCET-I contribution is given by the diagrams in Fig(2] Obviously, 
these contributions are leading order with respect to the strong coupling a s but subleading with respect 
to the electromagnetic constant a. Suppose that the incoming and outgoing quarks are hard-collinear. 
Then in the hard region where all momenta of the particles in the box diagrams are hard q^ 2 ~ Q 
this process can be described by perturbative QCD and we obtain the leading order contribution into 
coefficient function in front of certain SCET-I operators. But this region is only a part of the whole 
answer. For instance, one can also expect that the contribution from the region where only one photon 
is hard but the second is soft, for instance, q± ~ A, q% ~ Q can also be relevant. Therefore the SCET 
factorization must properly take into account all possible dominant regions associated with factorization 
of the hard modes. 

In the following we use two light-cone vectors n = (1, 0, 0, —1) and n = (1, 0, 0, 1) and choose the Breit 
frame where the initial and final proton are collinear to the z-axis. Then the light-cone expansion of the 
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external momenta read 

a = Q\-Q\, V^Q\, V'^Q\, (25) 

k^-Q^+ l -Q^ + k u V~± Q l + *Q?t + k x , (26) 
z 2 z 2 z 2 z 2 

where we introduced convenient dimensionless parameter z = —t/s and notation z = 1 — z. The Breit 
system is convenient for description of large Q 2 behavior of the nucleon FFs and it is natural to use it 
for the consideration of the TPE contributions too. 

Let us briefly describe the SCET notations used in our paper. We will use the SCET formulation 



in coordinate space developed in 43 . For the SCET fields we use following notations. The fields £ n , 



and denote hard-collinear quark and gluon fields associated with momentum p' and p, 



respectively, see Eq.(25). As usually, the hard-collinear quark fields satisfy 

iCn = o, - 0. (27) 

The fields q and Aft denote the soft quarks and gluons which also enters in the SCET Lagrangian. These 
fields describe particles with the soft momenta fc M ~ We also use the convenient notation for the 
gauge invariant combinations often appearing in SCET such as 

X „(Xn)= W n (An)Cn(An), (28) 

where the hard-collinear gluon Wilson line (WL) read: 



W n (z) = Pexp^igJ ds n- A {n) {z + sn)\ . (29) 

In the QED sector we also split the fields according to the same SCET prescription. Let us notice 
that the pure QED corrections (electron vertex and selfenegies) can be computed exactly and for such 
calculations such representations is not required and as we noted earlier we are not going to consider 
such corrections. But in case of TPE corrections the hard and soft photons are involved in the different 
dynamics and SCET representation is useful. 

In the wide-angle kinematics we have four directions with the large energy flow therefore it is useful 
to introduce two more auxiliary light-cone vectors associated with the lepton momenta: k and k' 



2kf „ 2k'v 



Q Q 
Then arbitrary vector can be decomposed as 



{v-v) = 2. (30) 



V i = {V-v)— + (V-v) — +V T . (31) 

Let us denote by subscript T the transverse components Vr with respect to v, V basis in order to distinguish 
them from the transverse components in n, n basis which we denote as V±. Using vectors v, v we introduce 
the hard-collinear lepton fields C and ( v 

fcv = 0, Kv = 0, (32) 

and the soft photon field . This is a minimal set of the auxiliary fields in QED which we need for 
our consideration. We do not introduce the collinear photon fields because they are not required for 
description of TPE amplitudes. We always use the Feynman gauge for the gluon and photon fields. 

We start our consideration from the one photon exchange. Let us mention that we consider now only 
the soft spectator scattering contribution. The soft spectator scattering for nucleons FF has been studied 

2 In Ref. |43| these modes are introduced as ultra-soft. In this work we use different terminology suggested in Refs. |45|46| 
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in 31 37 . Now we briefly repeat the results relevant for our calculations. To the leading order the SCET 
expression for the amplitude A ep can be written as 

Al v ~ ^ u(k')^u{k) C + (Q, M ) (p'|C^(0) b) SCET (33) 

where the leading power SCET-I operator reads 

01(0) = ^ e 9 {Xn{®)ll Xfi(O) + Xn(0hlXn(0)} , (34) 

9 



and denotes the hard coefficient function. Let us note that the operator in Eq.(34) has the additional 
term comparing to the operator considered in 31 . This term given by XnlxXn in Eq.(34| can be 



associated with the scattering of the hard photon on the antiquark. Such contribution describes the 
hadronization of the hard-collinear antiquark into proton. Because the proton is dominated by three 
quark state such term is suppressed by the powers of the hard-collinear scale in the transition from SCET- 
I to SCET-II. If the hard-collinear scale is not large and we restrict our consideration only by the SCET-I 
we can not neglect such operators. This moment has not been considered in our first paper [31]. The 



relative sign between the two terms in Eq.(34) is dictated by the C-parity. The subscript SCET in Eq.(33) 



(n) (n) 

means that the matrix element has to be computed with SCET Lagrangian £ SCET = £ SCET + £s CET + £ so [ t 
which represent by the sum of Lagrangians describing each collinear sector. The explicit expressions for 

The SCET matrix element 



£soet with the fields defined in position space can be found in Refs. 



43 



44 



in Ec[.(33) describes the interaction of the jets of the hard-coHinesT particles represented by Xn,n 

with 

the soft background described by the soft quark and gluons in the SCET Lagrangians. This dynamics 
can be also understood as a soft-overlap of the initial and final hadronic states. The parametrization of 
this matrix element reads 

(P'\ O" + (0) \p) SCET = N(p')^N(p) fi(Q,[i), (35) 

where /i is the SCET FF. By construction it depends on the hard-collinear scale therefore we indicate 
this writing in the argument Q, the second argument /i denotes the renormalization scale. Performing 
the matching it is convenient to put /i = Q and then to evolve it down to the values of the hard-collinear 



scale /i ~ ^AQ. The evolution is described by the renormalization of the SCET operator, see e.g. 47 
In what follow we suggest to use simple notation 

f 1 (Q lf i = Q)= fi(Q). (36) 

Taking into account that the tree level coefficient function reads 

C + (Q,n = Q) = l + G(a s ), (37) 

one obtains 



A % = u(k')ru(k) N(j/)^N(p) fi(Q). (38) 



Aira 

— „ vw , _ vv 4 
Comparing this with the A 7 in Eq.([7]) one finds (see details in Appendix A) 

G M (Q 2 ) = fx{Q), F 2 (Q 2 ) = — h{Q). (39) 

FF F 2 is suppressed by factor m 2 /Q 2 due to helicity flip. Moreover, this expression is not complete 
because it also includes the contribution from the subleading operator, see e.g. [3l], which we do not 
consider now for simplicity]^] 

3 In [3l| the kinematical power correction with SCET FF f\ was missed. The complete leading power contribution is 
given by the sum F2 = 4m 2 /i/Q 2 + m?CB * f'z/Q 2 where Cb is the subleading coefficient functions in notations [311. 
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The generalization of this approach to next-to-leading order in QED can be done following the same 
line. It is convenient to include the soft photon field into SCET Lagrangian similar to the soft gluon field 
and perform the expansion at the last step. Then for the elastic amplitude we can write 

A ep ~ (p\ k'\ T{C» V * 0%Ol + CT * OtOfi \k,p) SCET , (40) 

where the hard coefficient functions C± describe the hard subprocesses, the asterisk denotes the convo- 
lution integrals in position space and the SCET operators 0±" describe the dynamics of hard-collinear 
and soft particles and read: 

0£ = 0£(Ai, A 2 ) = Yj iXn(Xinhlxn(X2n) ± x s (A 2 n) 7 ^(Ain)} , (41) 

flavors 

O v e = 0»( m , m ) = Xv{vie)l u Xv{me), (42) 

where we introduced 

Xv =CvW v , X v = Wkv, (43) 

with the hard-collinear photon WLs 

Wy(z) =Pexp|-i e y dt v- B$(z + vt)\ , W v (z) = Pexp {~ ie J dt v - B$(z + vt)\ . (44) 
Taking into account that these WLs are not required for the TPE calculation we can skip them assuming 

oiimM-UmeWUme). (45) 

The similar simplification we also assume for the operators 0±. This implies that we reduce our 
consideration only to TPE diagrams in Figj2"P|The operators 0± represent the operators which arise at 



leading power and leading order in a s . The other operators will provide only subleading contributions 
according to SCET-I power counting or subleading in a s . 

The C-even operator O^L can appear only due to the hard two-photon exchange. Corresponding 
hard coefficient function CO" can be associated with the hard region in the TPE diagrams in FigJiJ i.e. 

~ 0(a 2 ). Hence computing the matrix element (O^Og) we can neglect by the soft photons in the 
SCET Lagrangian because they provide only higher order 0(a 3 ) contributions. Therefore the leptonic 
matrix element can be factorized 

(p', k'\ T{Ct v * 0»_O v e } \k,p) SCET = (k'\ 0» |fc) SCET * Ct v * (p'\ Ot b) SCET (46) 

The hadronic matrix element gives the new SCET FF which we define as following: 

(p'|0^(A l5 A 2 ) b) SCET = e *(f'^x-i( P n)A 2 N(p')tt^N{ V ) (47) 



The dependence on Q and fj, has to be understood in the same way as in case FF /i defined in Eq. ( 35 1 



Substituting this into Eq.(46l yields 



(p',k'\T{C» v *0>tO»} IMscet = ^C^(z,Q 2 ) u vl »u, N(p')<^rtN(p) 7i(Q,/i), (48) 

where C_ is the Fourier transformation of the coefficient function Ct!_ v defined in position space, u v ^ 
denote large components of the electron spinors = 0, pUy — 0. Notice that such approximations 
is exact in the limit of massless electron m e = 0. The factorization scale fip differs from the QCD 
factorization scale fi because it describe the factorization of the QED loops. In Figfl] the contribution 



given in Eq.(48l corresponds to the first graph on rhs. 



4 Let us note that each individual box diagram has the regions associated with collinear photons but in the sum these 
contributions cancel due to gauge invariance. 
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Consider now the first term in Eq. ( 40 ) . This contribution must be computed taking into account the 
presence of the soft photons in SCET Lagrangians. The situation in the leptonic sector can be simplified 
if we factorize the soft photons using the fields redefinition, see e.g. [48] 



(49) 



(50) 



where the soft photon WLs reads 

yj(0) =Pexp^-ieJ dt v ■ B (s) (to)j , 5 C (0) = Pexp dtv-B {s \tv) 

The new electron fields C^l do not interact with the soft photons anymore and therefore the leptonic 
matrix element is factorized 

(p>, k'\ T{C^ * 0%Ol} \k,p) SCET * Cf * (p'\ T{0» + {\ u A 2 )YJ (0)5,(0)} |p) SCET 

x (k'\ {e\mv)Yd 0) (v2v)} |fc) SCET , (51) 

and can be easily computed^ 



(k'UCKmvh^rimv)} I^)scet = ji^-to*™ u vl »u,. 



(52) 



The hadronic matrix element can be considered as a generalization of the leading order matrix element 
defined in Eq.(35) in the presence of the soft photons created by lepton source. Therefore we can write 



(p'\T{Ol(\ u \ 2 )Yj (0)5,(0)} \p) 



SCET 



LL(T 

91 



a 



fi(Q) + -9i(*,Q) 



+ vlv^_ aA ^^g 3 (z,Q)}. |.-,T, 



where the SCET amplitudes <7x,3 represent the next-to-leading QED correction arising from the interaction 
of the soft photons with the hard-collinear and soft spectator quarks. The coefficient 4 ^ 2 - 1 is introduced 
for convenience. These functions also depend on the factorization scales /i which we ignore in Eq. ( 53 1 for 
simplicity. Substituting (52 1 and (53) into (51) we obtain 



(j/, k'\T{Cl v * O^O"} IMscet = C£Wue Nip'rfllNip) 



LLC 

91 



h{Q) + -9i{z,Q) 



a4(z - 2) 



■v x vi — 



9?,{z,Q) 



(54) 



In Fig[T]the contributions with 51,3 in Eq.(54| corresponds to the second graph on rhs. The coefficient 
function is fixed by the one photon approximation. Using Eqs.(33) one obtains 



<T = !r-QrC + (Q,ri = gT^f + 0(a s ) + 0{a 2 ). 



(55) 



The terms of order 0(a 2 ) includes pure QED corrections like electron and photon self energies, photon- 
electron vertex corrections and also hadron vertex and self energy which we will not recompute in this 
work. Using Eq.(55), the TPE contribution to the amplitude AJ^ can be written as 



4na 

~Q 2 



u vl v u,N(p')L± 1 lN(p)\g^ fi(Q) + - gi (z,Q) + -K^P" -g 3 (z,Q) 



+ u vl »u, N(p')^ ± N(p) Ct v {z,Q 2 ) JKQ), 



(56) 



3 Here we do not describe the formal details related to transition from SCET-I to SCET-II in the leptonic sector. 
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where we used that 



«4^_2) = l^p,^ m ^^^(p). (57) 

" 7T z l 4 s 4 ~~ 



The expression in ( 56 ) is our final expression for the soft spectator scattering contribution describing 
elastic ep-scattering amplitudes with the TPE corrections as in Figj2j In order to use this formula in 
the phenomenological applications we have to compute the hard coefficient function and estimate 



the SCET FFs which enter in Eq.(56). 



3 Calculation of the hard coefficient function 

The calculation of the hard coefficient function can be done by matching of the box diagram in 
FigjZj In order to perform the matching we as usually consider quarks instead of proton states^] The 




Figure 4: Two-photon exchange diagrams with external momenta required for calculation of the hard 
coefficient function 



expressions for the diagrams with massless quarks are well known, see e.g. 38 49 and reads 



A% c = ^u vl ^ q n j 7 ^e 2 h + °^e\ / 8 1 q fi + v n {. . . }v n , (58) 

where and q n fi and i>„ jft denote the SCET spinors for the quark and antiquark, respectively. The ampli- 
tudes fi read (recalling z = —t/s) 

„ - , , s 1 , 2 1 z In 2 z z In 2 f 1 , ztt 2 n 2 

Re fx = nz n— + - n 2 z ^ hxz 1 , 59) 

J1 A 2 2 41-z 4 2 4 2' 17 

2 — z , n 2 — z , o z 1, , z 2 — z 9 /N 

Re/ 3 = — - In 2 z - -— In 2 - + - In z + In — ir 2 , 60 

2z z 2 z z z 2 

where A 2 denotes the artificial photon mass which plays the role of IR scale. To avoid misunderstanding 
let us note that A 2 is IR QCD regulator and is not related with the IR QED regularization. We also do 
not consider the imaginary parts because all discussed observables are sensitive only to the real part of 
the TPE amplitudes. 

On the other hand from discussion in the previous section it follows that corresponding amplitude 
can be written as 



A X = -q^u v Yu, ( q nl lq n + Vtftfc) l-g* (z, Q) + ~P v Klgl{z, Q) + C™(z, Q 2 ) Ff{Q)\ , (61) 



where FFs JF\ and g\ 3 denote the SCET amplitudes computed with quarks. FF T\ is easily computed 
from the tree level diagrams and reads 

T\ = e 2 + (D(a s ). (62) 

Calculation of g\ involves one-loop diagrams shown in Figj5j These diagrams must be computed in 
SCET-I which describes the coupling of the soft photon with the hard-collinear quarks. The one-loop 

6 More explicitly, we must consider quark plus antiquark state. However antiquark contribution can be easily restored 
from the quark calculation. We also consider one flavor with electric charge e q . 
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V . 




Figure 5: The one-loop diagrams a, b contribute to the SCET perturbative function g\ associated with 
TPE contribution. The dashed line with label v and v denote the soft WLs, the crossed quark lines show 
the different attachments of the soft gluon. The diagrams c and d are related to the vertex corrections 
and therefore were neglected. 



diagrams in Figj5]have UV-divergency which is regularized using dimensional regularization D = 4 — 2e. 
We present the details of this calculation in Appendix B. After UV renormalization the results reads 

gl(z,Q) = lnzln^f 7?, $£(*,Q) = 0. (63) 

At formal level such contribution arises form the mixing between the operators YJ S v O+ and Ot_ due to 
soft photon interactions. 

It is convenient to define the two independent coefficient functions Cm,3 as 

Ct"(z, Q 2 ) = gT-C M (z, Q\ p F ) + -P^K^-C 3 (z). (64) 



Using the explicit results for T\ and gf 3 and comparing two expressions for A^ 8 given in Eqs.(58) and 
(61 ) we obtain 



g f r-C M (z ) Q 2 , f i F ) + -P v K» x -C 3 {z,Q z ) = <?7 r /r-lnfln^- +-P»K£-f 3i (65) 

TT S TT TT \ X Z J S TT 

which yields 

1 z In 2 z z In 2 z/z 1 



2 41-2 4 2 42' 



C M (z,Q ,np)=]nzhij^- + -1n z -lnz — + (66) 

Cs(z, Q 2 ) = In 2 z - In 2 I + -_ lnz + In * - 2 —^tt 2 . (67) 

2z z 2 z z z 2 



Notice that the soft scale A 2 cancel in the final expression for Cm as it should be. In Eq.(66) we write 
ln[s//Xjr] with the total energy s for convenience. 

Therefore our final result for the soft spectator contribution with TPE amplitude can be written as 



1 a 



+ u v fuv N(p'y-f$ ± N(p) ~{g 3 (z,Q) + C 3 (z,Q 2 )F 1 (Q)}. (68) 

Let us now discuss the complete factorization formula describing the TPE contribution. We suggest 
that the TPE amplitudes in Eq.Q in the region of wide-angle scattering (23 1 can be described by 
following tentative factorization formula 

*3* a (e,Q 2 ) - 5G${e,Q 2 )+5G M \e,Q 2 ), (69) 
F 3 (e,Q 2 ) - F^(s,Q 2 )+F^ h \s,Q 2 ), (70) 
6F? c (e,Q 2 ) = SF^(e,Q 2 ) + 6Ft\s,Q 2 ), (71) 

where the indices (s) and (h) denote the contributions related to the soft and hard spectator scattering, 
respectively. 



12 



The hard spectator contributions describe the scattering where the all spectator quarks are involved 
into the hard subprocess. Schematically these contributions are shown by the third diagram on rhs in 
Fig{T] Schematically one can write 



ff } = *( Xi )*H 3 (z,Q 



5G$ = 9( Xi )*H 



M 



(z,Q 2 ; XuDi) * 



(72) 
(73) 



puted in Refs. [32 33 
Fig{2j and they are of order a s 



where symbols H 3 m denote the hard coefficient functions and the asterisk is used as notation for the 
convolution integrals over the collinear fractions Xi,yi. The leading order coefficient functions H 3 m com- 
The hard spectator contributions are dominated by one hard gluon exchange, see 
The nucleon distribution amplitude ^f(xi) describes the nonperturbative 
overlap of the collinear quarks with the nucleon state. The hard spectator contribution to the helicity 
flip amplitude i*2 can & l so be written in the similar way but it is power suppressed and the convolution 
integrals ill defined due to the end-point singularities. Therefore we will not consider this FF in the 
present publication. 

The soft spectator contributions describe the scattering where the spectator quarks are soft. These 



TPE contributions have been discussed above. Comparing Eq.([8|) with expression in Eq.(68) (the match- 
ing is quite similar to one described in Appendix A) we obtain 



5G$(e,Q 2 ) = -{ g 1 {z,Q^ F ) + C M {z,Q 2 ^ F )F l {Q)} 



v p (s) _ a v 

o 3 — 
m z TT s 



{g 3 (z 1 Q) + C 3 (z,Q 2 )F 1 (Q)} 



(74) 
(75) 



where we assume z = z(e) and v defined in Q. 

The soft contribution to the helicity flip amplitude is also suppressed by the power of large scale l/Q 2 
and can be written as (for simplicity we do not write the arguments on rhs) 



SF^(e,Q 2 



8F, 



a 4m 



2 r 



abl TT 



.9i + Cm J- 1 + -{C 3 F 1 +g 3 } 
s 



(76) 



In this expression we computed only the kinematical corrections similar to Pauli FF F% in Eq.(39). The 



complete answer includes also contribution with the matrix elements of the subleading SCET operators 



SR 



(<0 



subl 



which we do not consider for simplicity. 



As we can see from Eqs.( 74|75 ) the two leading power amplitudes are described by the three nonper- 
turbative functions: SCET FF J-\(Q) and SCET amplitudes 31,3(2, Q)- If both photons are hard then the 
soft spectator scattering is described by the pure QCD sector. In this case the hard-collinear dynamics 
involves only the one hard scale Q and described by the SCET FF F\(Q). Therefore the e-dependence 
is described by the hard coefficient functions C 3 _m and can be computed in perturbative QCD. The 
situation is changing if the one photon is soft. Such soft photon interacts with the hard-collinear and 
soft constituents and therefore the soft-overlap contribution is more complicate and described already by 
the SCET amplitudes 171,3(2, Q 2 )- Remind that we assume that the hard-collinear scale is not large and 
we can not use it for perturbative calculations. In this case the e(z)-dependence originates from the soft 
dynamics and can not be computed from pQCD. 

Another difficulty is related to the separation of the amplitudes on the soft and hard spectator 
contributions. The simple sums as in Eqs.( 69fl7lj are motivated by the structure of the SCET-I operators. 
However the overlap of the soft and collinear sectors in SCET-II makes such separation ambiguous. As a 
rule this leads to the end-point singularities in the collinear integrals defining the hard spectator terms. 
In Ref. 



37 



it was discussed that such situation is actual even for the FF Fi where one could expect that 
the hard spectator contribution is well defined. In general, the separation on the soft and hard spectators 
contributions can be well formulated only within a certain regularization scheme which allows to treat 
the soft and collinear sectors separately and consistently. From this point of view our results are still not 
complete because we did not provide such a separation scheme. This question we will discuss in the next 
two sections. 

Finally let us note that the QED dynamics considered here is closely connected with the underlying 
QCD dynamics. In particular the energy and virtuality of the soft photon can not be smaller then the 
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typical QCD soft scale A, because otherwise the photon interacts with the total proton charge and can 
not resolve the constituent charges inside the target. Therefore our QCD IR regulator A introduced for 
convenience as photon mass must cancel in the hard coefficient function which describes the scattering 
of the hard photons. The standard QED IR-divergence is included in the SCET amplitude <?i which is 



defined via the matrix element in Eq.(53). 



Let us also note that each separate diagram in Fig|4]has also the contributions associated with the 
collinear regions (one of the photons is collinear to one of the external momenta) . Moreover these collinear 
regions overlap with the soft one. However in the sum of the diagrams, the collinear contributions cancel 
and this cancellation is required by gauge invariance because the electric charges of the quark and lepton 
are different. For the interested reader, we discuss this in more detail in Appendix C. 

In conclusion, we suggested the factorization formula for the TPE soft spectator scattering contribu- 
tion. We performed the matching and computed the coefficient functions in leading order in approxi- 
mation a s . In order to apply this result to the phenomenological analysis, we need to define the SCET 
amplitudes which arise in our approach. 



4 Estimate of the soft spectator contribution with T\ from wide- 
angle Compton scattering data 

4.1 QCD factorization for wide-angle Compton scattering process 



The SCET FF J-\{Q) defined in Eq.(47) describes the long distance interactions in the soft spectator 
contribution with the hard photons. In order to estimate this quantity at the different values of Q 2 
one can use the universality of its definition. In other words if one can find the other process which is 
described by the same SCET FF then one can use the corresponding observables in order to estimate 
J-\(Q). In our case such analysis can be carried out for the wide-angle Compton scattering (WACS). The 
cross section of this reaction and polarization asymmetries in the region of large Mandelstam variables 
(large scattering angles 8 C 



90° in the cms frame ) have been measured in JLab 50 



The theoretical description of this reaction has been presented also long time ago 52 53 and it is 



based on the dominance of the hard spectator scattering which leads to the collinear factorization for 
the dominant amplitudes. In this case the situation is very similar to the situation with the nucleon 
electromagnetic FFs. At large energy and momentum transfer s ~ — t ~ — u 3> A 2 the leading power 
contribution to the WACS amplitude Tj can be described as collinear convolution of the hard coefficient 
function Hi(s,t) with the nucleon DAs \& 



if t) = * H^s, t- x i>yi ) * ¥( Xi ) 



(77) 



where the asterisk denotes the convolution integrals with respect to the collinear fractions Xi,yi. The 
reduced diagram describing the hard spectator scattering is shown in Fig J6) But at the same time one 



Figure 6: Reduced diagram describing the hard spectator scattering contribution in the WACS. 
can easily find that the soft spectator scattering can also contribute to the same power and that this 



situation is quite similar to the one discussed for the FF F\ in 31 351. For instance, one can investigate 



case of FF F x in 31[ 35 
momenta k 



the soft region for the diagram in FigjT] in the same way as it was done for the analogous diagram in 

Such analysis allows to see that the contribution from the domain where 
A are soft is suppressed by the same power of Q 2 as from the domain where 
2fi ~ ~ Q- This observation provides the arguments that the complete 
factorization is described by the soft and hard spectator contributions where the last one is described 



>2/, 



these momenta are hard: k 
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Figure 7: Example of the diagram which has the leading power contribution from the region with the 
soft spectators fc 2At ~ k 3 ^ ~ A 



by the collinear factorization suggested in [521 . Moreover phenomenological analysis of the experimental 



data 50 51 allows one to conclude that the dominant contribution in the cross section is provided by 



the soft spectator scattering mechanism. The best description of existing WACS data has been achieved 



using the so-called handbag or GPDs approach 54-57 which can be considered as model for the soft 



spectator scattering. However, the systematic analysis within this approach is rather problematic because 
the parametrization of the soft overlap contribution by the GPD matrix element is not consistent with 
the IR-structure of the QCD diagrams, see for instance [56] , 

In the current section we consider the description of the WACS in the SCET framework. Such 





Figure 8: Schematical representation of the WACS factorization. The soft-overlap contributions are 
described by the SCET FFs T 1 and Q x . 



consideration is quite similar to the analysis of the nucleon FFs carried out in 31 . We expect that at 
the region s ~ — t ~ — u ^> A 2 the WACS amplitudes T can be described within SECT-I approach by 
the following tentative factorization formula 



T i (s > t)=T^(s > t)+T^(s ) t), 



where the hard spectator contribution T^ h > is given by Eq.(77) and the soft spectator term read 
if } ( s , t) = Cf(s, t, /u)7i(<2, fi) + C?(s, t, /j,)Gi(Q, m), 



(78) 



(79) 



where the hard coefficient functions C\ ,a describe a hard subprocess, the functions T\ and Q\ describe 
the hard-collinear and soft interactions and can be defined in the framework of SCET-I approach. Wc 
provide their explicit definitions below. In Fig(8]we illustrate the WACS factorization in terms of reduced 
diagrams. 

In order to proceed further let us introduce following definitions. The amplitude of the process 
j(k) + p(p) —> 7(A/) + p(p') can be written as [58] 



(p', k' out\ in k,p)=i (2tt) 4 S(p + k + p + k') e*" (k r ) e" (fc) e 2 T" u , 



with 



TV = i J d A x (p'\ T{J£ m {Q)r em {x)} \p) , 

where J£ m represents the quark electromagnetic current 

=J2 e 1 ^(x^qix). 



(80) 



(81) 



(82) 
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The hadronic matrix element T^ v is parametrized in terms of six independent scalar amplitudes 7$ as 



p/fl plV AT /J ATI/ 

=N(p') i _ (T X + ^ T 2 ) - — _ (T 3 + £ T 4 ) 



p/2 



P' 2 7C 2 



iV 2 

P' 2 X 2 



z 75 |T T 6 I AT(p) 



(83) 



In the last formula we used 



N» = S [»a0j}P a ^(p-p'fK'r, e i23 = +1, 



s= (p + fc) 2 , t={p' -pf 



— ~ -(1 - cos6» cm ). 
s 2 



(84) 



One can see from Eq.(83) that three amplitudes 72,4.6 describe the Compton scattering without nucleon 



helicity flip and the three other with helicity flip. 

We again consider a Breit-like frame and therefore the light-cone expansions of the momenta can 
be easily obtained from Eq. ( 26 ) substituting the lepton momenta k and k' by photon momenta q and 
q', respectively. Therefore the formulation of the appropriate SCET degrees of freedom in the collinear 
sectors associated with n and n light-cone directions is the same as discussed in Section [2j 

In the following we only discuss the factorization for the soft spectator contribution. The SCET-I 



factorization implies that the T-product of the electromagnetic currents in Eq.(81 ) can be represented as 



the convolution of the coefficient function with a certain SCET-I operator constructed from the SCET-I 



fields £ n , 



00 



(n) 



This schematically can be written as 



T{J? m (0)J» m (x)} 



(85) 



where the asterisk denotes the convolution integrals in the coordinate space. We assume that the dom- 
inant contribution in the wide-angle kinematics is provided by appropriate operators with the minimal 
dimension according to SCET-I counting rules, see Refs. 39-44 . Further in our analysis of the soft 



spectator scattering contribution we consider as a leading power terms only the contributions associated 
with the SCET-I counting, i.e. we consider the powers of the large scale Q arising from factorization of 
the hard modes. We do not consider the full power counting analysis which includes also the SCET-II 
counting. We assume that our consideration is actual for the region where the hard-collinear scale is still 
small and therefore we can not perform expansion with respect to this scale. 

In our case we find that the set of the suitable operators which is consistent with the symmetries and 
has the minimal SCET-I dimension consist from the quark operator 0"(Ai,A 2 ) defined in Eq.(41| and 

from pure gluon operator 0"^(Ai, Aa) which includes the transverse gluon fields A^]_, 



written as 



A^_ and can be 



'(Ai,A a ) = Tr{[iyt(A 1 n)/9 ct± I¥„(A 1 n)][^(A 2 n)^ ± ^(A 2 n)]} 



(86) 



where the covariant derivative iD^ = id^ 



gj^ji^n) ac ^ g Qn fjgidg on iy inside the square brackets 



Below we restrict our discussion only by the leading order in a s . Then we can skip discussion the gluon 
operator because the corresponding hard coefficient function arises only at the next-to-leading order. 
Therefore at leading order we have contribution only due to the quark operator O" (Ai, A 2 ) that yields 



T{JL(0)Jem(x)} * Cr * O^Ax, A 2 ). 



(87) 



Substituting this into matrix element in Eq.(|81[) and using Eq.(47) we obtain 
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+ crossed 



Figure 9: Leading order diagrams to the hard coefficient function describing the soft spectator contribu- 
tion in WACS. 



where C^ va (s,t) corresponds to the coefficient function in momentum space. From Eq.(88) we can 



conclude that all six amplitudes Tj introduced in Eq.(81) are described in term of one SCET FF T\. 
However helicity flip structures in this case can be obtained only as subleading contributions. We also 
observe dependence on the total energy s is completely described by the hard coefficient function. 

In order to compute Cfi va one has to perform the matching for the tree level diagrams shown in Figj9j 
We skip the details of this relatively simple calculation and provide the results for the WACS amplitudes 



defined in Eq. ( 83 ) . For the scalar amplitudes we obtain {Q 



-t, m is nucleon mass) 



W _ _ *(2s + t) 

2 4 ~ *(* + *) 



w _ 



s(s + t) 



(89) 



r ( 8 ) = m 1 (2s 



tf 



t2 s(s + t) 



(*) 



W = ml {2s + tf 
3 i2s(s+t) 



Fi{Q) + % 



From these equations one can conclude that helicity flip amplitudes T[ 3 5 are suppressed by power of 1/Q 
(to see this explicitly one has to pass to dimensionless amplitudes in Eq.(83)). The terms Tj s ' appearing 



(s) 



(90) 



on rhs of Eqs.(90) denote contributions associated with the subleading SCET operators which has not 



been considered in our analysis. 



The separation of the hard and soft contributions as described in Eq.(78) is not simple because it 



usually implies an additional regularization. Such regularization is required in order to separate the soft 
and collinear sectors and to describe correctly the so-called large rapidity logarithms, see for instance 
Refs. 59-62] and references there in. Such problem is clearly seen when one computes the hard spectator 
contributions and obtain that the collinear convolution integrals are divergent in the end-point region. 
This is a clear indication that one has to consider the overlap of the soft and hard spectator contributions. 

In Ref. [37] the soft and hard spectator contributions describing the FF F\ has been investigated and 
it is demonstrated that in this case the problem of the soft-collinear overlap also exists. We suppose 
that collinear factorization describing the hard spectator contribution in the TPE and WACS case is 
also violated due to specific end-point divergencies. Therefore we suppose that the most optimal way to 
proceed further is to explore the so-called physical subtraction scheme suggested in 63|64 for description 
of different £?-decays processes. 

The idea of this approach is to exclude the FF T\ from the expression for the physical amplitudes 



and write the relations between the physical amplitudes directly. In our case, using Eqs.(79) (neglecting 



the gluon contribution at leading order) and Eqs.(77) we can write 



(91) 



C 2 {s,t) VtfV C 2 {s,t) 
Such formal expression implies that we use some regularization in order to define the divergent quantities 



in both sides of (91). Substituting expression (91) into Eqs.(79l for T^ q we obtain 



Ti{s,t) 



C i (s,t)^^- + ^(y i )*{H i (s : t ] Xi ,Vi) 
C 2 {s,t) 



H 2 {s,t; 



> Vi 



C 2 (s,t) 



9(xi), £ = 4,6. 



(92) 



The end-point singularities arising in the collinear convolution integrals \& * Hi * \& must cancel in the 



expression in Eq.(92) because the other terms with the physical amplitudes do not have any problems. 



Hence we can study the relations between the physical amplitudes T{ systematically order by order in 
a s . Taking into account that the coefficient functions Hi are of order a 2 s and therefore assuming that 



17 



corresponding corrections in ( 92 1 are small we obtain the leading order relations 



r 4 (8,t) 

T e (s,t) 



C 4 (s,t)%^ = -T 2 (a,t) + 0(a a ), 
C 2 (s,t) 



C 6 (s.t) 



T 2 (s,t) = t 
C 2 (s,t) 2s + t 



T 2 (s,t) + 0(a s ). 



(93) 
(94) 



where the coefficient functions 6*2,4,6 have been computed above to the leading order accuracy. Let us 
also note that the next-to-lcading corrections to these relations arise from the corrections to the 6*2,4,6- 
The nice feature is that SCET FF T\ does not depend on s. Therefore this allows one to substitute 



Eq.(91) into other expressions with T\ even for the different values of s. Such situation is important in 



our case because the cross section for the elastic ep-scattering have been measured for different values of 
s. 

From Eqs.( 93|94 | one can easily find that to leading order accuracy the ratios Tj/Cj satisfy 



T 2 (s,t) r 4 (a,t) T 6 (s,t) 
C 2 (s,t) ~ Ci{s,t) ~ C 6 (s,t)' 

Therefore it is convenient to introduce following notation 



1l{s,t) 



T 2 {s,t) _ 2s + t 
C 2 (s,t) ~ s(s + t) 



T 2 (s,t) + <D(a s 



(95) 



(96) 



Contrary to the SCET FF T\ the ratio TZ is a well defined object free from the rapidity regularization 
dependence^ This ratio formally depends on the two variables s and t but according to construction, in 
the region where QCD factorization is applicable and power and higher order in a s corrections are small 
the dependence on s must be small, i.e. we must observe that 



dK(s,t) 
ds 



0, 



(97) 



which is a good criteria of the applicability of our approximations. The ratio TZ for the different values 



of s can be extracted from the cross section for WACS measured in JLab 50 . Then this ratio can be 



used for the calculation of the TPE amplitude defined in Eqs.(74 75) 



4.2 Phenomenological analysis of the WACS observables and extraction of 
the ratio TZ 



The unpolarized cross section describing Compton scattering reads 58 

da ira 2 
dt (s — m 2 ) 



;W 00l 



(98) 



with 



Woo = g ( s ~ ™ 2 )(™ 2 - «)(l r 2| 2 + |T 4 | 2 ) + (m 4 - S u)|T 6 | 2 

+ m(s - u) Re [T.T; + T 3 T 4 *] + ^(4m 2 - t)(l^i| 2 + |T 3 | 2 ) - t\T 5 \< 



(99) 



Substituting the expressions for Ti into ( 99 ) and neglecting by the subleading power corrections we obtain 

2ira 2 



da 

'dt 



(s — m 2 ) 2 \ 1 — z 



1 



da 



KN 



dt 



(100) 



where da™ is the Klein-Nishina cross section on the point-like massless particle. This formula is very 
similar to one obtained in the handbag approach |54| except for the definition of the form factor. 



We assume that the coefficient function C2 in Eq.l96k is computed at factorization scale fi 
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There are two more observables which describe the correlations of the recoil polarization with the 
polarization of the photons 



K LS = 



(101) 



where ||(_L) refers to a longitudinally (transversely) polarized nucleon and R(L) denotes a right (left) 
handed photon. For the first one we obtain (again neglecting the power corrections) 



K 



LL 



2z(2-z) 
(2-z) 2 + z 2 



(102) 



In this case the SCET FF cancel and we simply reproduce the Klein-Nishina for the point-like massless 
particle. The second asymmetry Kls is sensitive also to the nucleon helicity flip amplitudes T\ 3 5, for 



which we used the approximations with the kinematical contributions only, see Eqs.(90 1. Then the leading 
order expression reads 

z~ (2 - z) (3z - 2) 



K 



LO 

LS 



1 



(2 - z) 2 + z 2 



where m denotes the nucleon mass. These asymmetries have been measured in JLab 
and t = -4.0 GeV 2 
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Kf% = 0.678 ± 0.083, 



Kfg = 0.114 ±0.078. 



Using our approximate formulas ( 102 ) and ( 103 ) we obtain for these kinematics 

0.70, Kl% 



0.066. 



(103) 
at s = 6.9 GeV 2 
(104) 

(105) 



The expression for the longitudinal polarization is in a good agreement while the transverse polarization 
is about factor two smaller but still within the error bars. Notice however, that in this kinematics the 
value of the u ~ —1.14 GeV 2 is not large comparing to the nucleon mass that potentially can provide the 
relatively large power corrections. Therefore the good agreement with the experiment numbers can be 
probably explained by lucky compensation of these dangerous corrections in higher orders. Nevertheless 
this analysis of the recoil polarizations allows us to make the following qualitative conclusion: it seems 



that the symmetry relations (93) and (94) provides the reasonable approximation for WACS. Therefore 



it allows us to use the data for the unpolarized cross section and formula ( 100 ) in order to extract 
information about the function 1Z. 

In order to extract the values \TZ\ from the data we take the data for the WACS from 50 j for which 
\t\ > 2.5 GeV 2 and \u\ > 2.5 GeV 2 . Then for these points we extract the values 



\K C 



lda cx P(s,t) 



(106) 



In Fig 11 we show the results for the |7?. exp | extracted from the data at three different values of energy: 
s = 6.8,8.9 and 10.9 GeV 2 and 2.5 < -t < 6.5 GeV 2 . As one can see from this plot obtained \K exp \ 
do not show any dependence on the values of s as required by factorization. Therefore we fit the Q 2 
dependence of the data by the simple empirical function 



n{Q) 



In 2 Q 2 /A 2 ' 



(107) 



which gives c = 0.10 ± 0.01 and A = 1.08 ± 0.03 GeV. 

In order to see the sensitivity of our approach to the possible power corrections we also perform the 
same extractions keeping the exact kinematics in the formula for the cross section. In this case we see 
that obtained |7\L cxp | are larger, especially in the region of smaller —t and the extracted values also show 
higher sensitivity on s. The solid line in this panel shows the same fit ( 107 1 which in this case yields 
c = 0.13 ± 0.01 and A = 1.09 ± 0.02 GeV. From this analysis we can conclude that the data at smaller 
values of \u\ or \t\ potentially may obtain sizable corrections from the power suppressed contributions. 
However because our calculation is incomplete this estimate can only be considered as a qualitative 
approximation. In our numerical calculation of TPE amplitudes we will use the results obtained without 
subleading power corrections. 
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• s=6.8GeV 2 
a s=8.9GeV 2 
■ s=10.9GeV 2 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 

-t, GeV 2 



Figure 10: The ratio 1Z extracted from the WACS data (50| and the fit (1071 without kinematical power 
corrections . The gray band shows the la error bands. 




Figure 11: The same as in Fig 10 but including the kinematical power corrections. 



5 The amplitudes in the hadronic model 

Unfortunately the extraction of the amplitudes gi t s from the other process at present is not possible. 
Therefore one has to build a model in order to estimate these quantities. 

As we discussed above, the splitting of physical amplitudes on the soft and hard spectator contribu- 
tion implies a certain rapidity regularization which allows one to consider the collinear and soft sectors 
independently. In the present paper we do not consider such scenario for g\ 3 and perform the numer- 
ical estimates of the hard spectator contributions defined in Eqs.( 73|72 ) using the nucleon distribution 
amplitude which does not provide the end-point singularities in the collinear convolution integrals. This 
allows us to consider the SCET matrix element which defines the amplitudes in (53 1 as a well defined 
quantity. 

Th SCET amplitude g\ depends on the factorization scale [ip. This factorization scale arises from the 
factorization of QED box graphs where [ip separates the the hard photon configuration from the regions 
with soft virtualities. The evolution equation with res pect to this scale can be obtained from the QED 
renormalization of the SCET operator YJS v O , f, see Eq.(63). Or equivalently the same can be easily fixed 



from the requirement that the physical amplitude does not depend on [ip . Using Eq. ( 68 1 yields 



fj, F 



-^—Al pc = ^ ^F-^—gi{z,Q^p) + ^F-^—C M (z,Q 2 ,li f )F 1 = 0. 
dfiF u^f dfip 



(108) 



From the expression for the Cm in Eq.(66) we obtain 



[If-, — gi(z,Q,fi F ) = 21nz Fi(Q). 
d\iF 



(109) 
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Solution of this evolution equation reads 



9i(z,Q,n F ) = gi(z,Q,Ho) + lnz In— |J"i(Q), (110) 

Mo 



where g\(z,Q, /j,q) denotes some initial condition. The logarithmic contribution in Eq.(llO) can be ab- 
sorbed into coefficient function Cm- Then in order to estimate the contribution associated with the 
amplitude g\ we have to fix fip = fj, and evaluate gi(z,Q, /Iq). 

Let us consider first appropriate choice of the scale hq. As we have mentioned before, factorizing 
the hard modes we pass to the SCET-I which describes the hard-collinear dynamics. Therefore one can 
expect that it is natural to assume that fig ~ y/ AQ. However the contributions associated with the 
hard-collinear photons cancel in the sum of TPE diagrams. This allows us to assume that factorization 
scale fj,p separates the hard and soft electromagnetic configurations in the box diagrams. Using this 
interpretation we can fix the scale to be of order A. With such choice <?i(/io — A) describes the TPE 
configurations where the soft photon carries momentum qi^ < A. Notice that in this case the virtuality 
of the involved active quark remains hard-collinear as it should be. Such choice for /ip explains its 
decoupling from the QCD factorization scale fi ~ y/AQ which enters in the definition of the QCD SCET 
FFs. Let us also mention that choosing /io — A one obtains the large logarithm in the coefficient function: 
Cm (a* f — f^o) ~ J^iln-Z However the combination — In is still quite small for the practically 

fj, r. M 

relevant values of s therefore we do not need to sum such contributions. 

Taking into account that the amplitude g\ describes the highly asymmetrical configuration where one 
of the photons carries the momentum q ifl < A one can conclude that it must include also the IR-sensitive 
QED contribution. Such term is very important for the correct treatment of the QED IR-divergencies. 
In what follow we suppose that the soft photon virtualities which are relevant for the QED IR-divergent 
part are much smaller then the QCD soft virtualities of order A: qf <C A 2 . Such ultrasoft photons 
can interact only with the point-like proton and their dynamics can be described in the framework of 
some effective theory with the nucleon degrees of freedom. These arguments suggest that we can try to 
estimate gi(no ~ A) using a model with a nucleon and excited resonance states. We may expect that 
such model can provide a reasonable estimate if the dominant contribution to 31,3(^0 ~ A) arises from 
the region of ultrasoft photon. From this point of view such calculation is quite similar to the hadronic 
model calculations with the one important difference: the virtuality of the soft photon is constrained by 
fj,Q ~ A 2 and the hard dynamics is factorized into the hard coefficient function C+. 

In order to perform our calculation we proceed as follows. First, we insert the intermediate nucleon 
states approximating the matrix element 



(p'\ T{Ot, (YjSv - 1)} | P ) SCET ^ f-^Y. V\ °+ ^ S >SCET {r, «l (tfS, " 1) \P) 



^ Y (P'\ MS, 1) |r, s) <r, S | Ol \ P ) SCET (111) 



where R denotes the nucleon resonance state and s describes its polarization. In Eq.(lll) we wrote the 
combination YJSv — 1 in order to stress that we consider only the next-to-leading QED contributions. 
Let us again to mention that we consider only the TPE interactions of the soft photons with hadrons 
neglecting the electron and hadron vertex contributions. The matrix elements with the soft photon WLs 
can be computed in the effective theory with the nucleon degrees of freedom. 

(r, a\ (YjS, - 1) \p) ~ J ' dx (r, s\ T { (yjff - l) , J%(x)teB^ (x)} \p) (112) 

where Jr(x) denotes the electromagnetic current of hadron R. Contracting the photon fields we obtain 

t{(YJSv-1) ,ieB^(x)} ~T^-iej db v • (to) - ie J dt v ■ B^(tv), ieSW(x)| (113) 

= e 2 fJ^eWJ-*L J !g + ^ ] (114) 
' (2tt) d [l 2 + ie] \[-{vl) + ie] ^[(vl)+ie] r 
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Substitution (114 1 in (112) gives 



(r,a\(Y^S v -l) \p)ce 2 (r,s\J%(0) \p) 



d D l 



- ie\ 



5 D {r-l+p) 



+ 



(vl) + ie] [(vl) + ie] 



(115) 



The similar expression ca n be obtained also for the second matrix element in Eq.(lll|. Substituting 
these expressions into Eq.(lll) and integrating over dPr yields 



(p'\t{o^(y:s,-i)}\ p ) scet 



e 2 u 2e 



d D l 



(2n) D J [I 2 + «e] I [-(vl) + ie] [(vl) + 



ie 



J2 m rV)> ( 116 ) 



where 



(p + l) -m 2 R 



(p + l,s\J%(0) \p) 



+ £<PV*(o) tf-ho) — 



{jpf -l) -m\ 



(p'-l,s\0>i \p) 



SCET • 



(117) 



This representation can be easily understood as calculation of the diagrams in Fig jl2| We use dimensional 
regularization with D = 4 — 2e in order to regularize UV-divergencies, see further discussion. 




P "> \_ P + 

J R R 




+ 




+ 




Figure 12: The diagrams in the low energy effective theory describing expression for the g±. The simple 
dashed line denotes WLs associated with vectors v and v, the black square denotes the vertex with 0+ 

In order to describe correctly the region associated with the small momentum I (soft photon) we need 
to perform in Eq.( 1171 some approximations. Namely, we neglect by small momentum I in the numerators 
and denominators in Eq.(117). 

(p'\ Ol \p + l,s)~ (p'\ Ol \p, S ) SCET , (p + l, s\ Jg(0) \p) * (p, s\ Jg(0) \ P ) . (118) 

i i i i , 

(119) 

(120) 



(p + l) 2 -m 2 R P+l- + m 2 - m 2 R (p> _ if _ m | -p'_l+ + m 2 - m\ .' 
In these formulas we do not skip the masses because 

p+l- % AQ, m 2 — m 2 R ~ (m — m R ) (m + m R ) ~ A (m + run) , 

and therefore the difference n 



— m 2 R is potentially comparable with the hard-collinear scale in the 
region of moderate values of Q where AQ ~ m 2 . Moreover this difference of the masses provides also a 
natural regulator for the QED IR-divergency. In this case only the elastic contribution has IR-divergency. 
Therefore we can write 



e 2 u 2e 



(p'\oI(y:s,-i)\ p ) scet - 

(2tt) 



d D l 



I 2 - \ 2 + ie] \ [-(vl) + ie] [{vl) + ie] 

x (M£<*(l)+M£*(l) + ...) 



(121) 



where we introduced the IR-regularization by the photon mass A. In what follow we restrict our considera- 
tion computing only the elastic contribution M£ a . This contribution is reqiered for the correct treatment 
of the QED divergencies in the physical cross section and usually it provides the largest numerical effect. 
Using (35) and 

(p\ J a (0)\p)~N(p) 7 a N(p), ^7V a (p, S )A> /3 (p,,s) = (p + m) Q0 , (122) 
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we obtain 



Mr 



N(p>) 



hi h{Q)}i^ F T + T 



p + l--\-ie 1 1 I — p'_/++ie 

Nip') ^Y ± N(p) h[Q) 



l hi h(Q) 



N(p) 



l--\-ie — l++ie 



Substituting this expression in Eq.(121| we obtain 



(p'\ 0* (YjS, - 1) b) SCET ~ e 2 N(p')^N(p) h{Q) J el , 



with the scalar integral 



d D l 



[l 2 -X 2 + ie] \ [-{vl) + 



l£\ 



[(vl) 



ie\ 



(In) 



■ is 



(In) 



■ IS 



(123) 
(124) 

(125) 

(126) 



This integral is the same as for the quarks in Appendix B. This is not surprising because we consider the 
soft photon exchange between the highly energetic particles (lepton and proton). The effective action 
describing such a process is again can be written in the the form of the soft photon WLs. The only 
difference is in the electric charges: the soft photon interacts with the total proton charge. Therefore for 
the real part we obtain 



Re J e , 



4tt 2 



In 



(vn) 



(vn) 



1 

^UV 



In 



A 2 



where u uv is the UV-renormalization scale in M S'-scheme. Obtained result can be rewritten as 



In 



(vn) 



(vn) 



In 



2(kp) 



2(kp> 



In 



In 



where we introduced s 



-u + in 2 . After UV-renormalization we find 



( P '\ Q» (YjS* - 1) \ P ) 



SCET = " N(p')^N(p) fl (Q)\n^\n 



(127) 



(128) 



(129) 



Let us discuss obtained result. The UV-divergency arises as soon as we neglect the small terms I 2 in 



the denominators of the nucleon in Eq.(119). Such approximation can be understood as transition to 
the effective theory which describes the interaction of the energetic particle^] with the soft photons. 
Appropriate effective action can be naturally described by the soft photon WLs directed along the proton 
momenta. Formally the UV-divergency is related to the multiplicative renormalization of WLs. The UV- 
renormalization induces the dependence on the scale u uv which we can be understood as highest possible 
virtuality of the soft photon consistent with the made approximations. It is naturally to expect that we 
can put y« uv = (j,q ~ A. In the final expression (129) we keep the exact expression In | -S- 1 because this 
expression provides the exact cancellation of the IR-divergent contribution. This moment is discussed in 
details in the next section. 

Our calculation has many common features with the soft approximation made by in the work [10] . 
In that work the TPE diagram with the nucleon intermediate state has also been computed neglecting 
the small photon momentum I in the numerator of the integrand (including the FF of the hard vertex) 
and in the hard photon propagator but not in the denominators of the electron and proton. As a result 
the integrand is reduced to the three-point function. Obtained integral is UV-finite but the answer 
unavoidably includes the region where the virtuality I 2 of the photon (and other particles) can be large 

I 2 ~ Q 2 . But such description of this hard region is incorrect because the dynamics of the hard particles 

I II 1 \ 2 
must be described using QCD. As a result the obtained in 10 38 result includes large logarithm In — 



\2 

£ tpo [Tsai] - In — In 
s 



G 



A I 



(130) 



'Remind that in our frame the leptons and protons move with the energies of order s ~ Q. 
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because in this case the hard scale s provides the required "UV cut-off" . Therefore we expect that such 
description of the TPE contribution in the region of large energies is incorrect in spite of the fact that it 
correctly treats the QED IR behavior. 

In our calculation we must not have any region with the hard virtualities because we already factorized 



the hard dynamics. Therefore it is important to perform the expansion in small momenta / as in (119) 
in order to avoid the double counting of the hard region. Then we do not obtain any large logarithms. 
The uncertainty in the choice of the renormalization scale fi uv ~ /io has a clear physical meaning and 
is associated with the applicability of the our soft approximation and also reflects an ambiguity in the 
calculation of the initial condition gi(z,fi Q ). 



Hence assuming the dominance of the elastic contribution and comparing Eq. ( 129 1 with ( 53 ) we obtain 



(X 2 

9i (z, Q, Mo) - - h (Q) In ^| In ^ 



- Gm{Q 2 ) In — 2 



g 3 {z,Q)~0. 



(131) 



where we used Eq.(39l to re-express SCET FF /1 in term of physical FF G 



We obtained that the second SCET amplitude g 3 in this approximation is very small. We suppose that 
this is the feature of our model in which the soft photon interacts with the hadron as with a structureless 
point-like particle. We expect that contributions with the higher spins in Eq.(121| also do not provide 



any significant contributions to g 3 . Probably, this amplitude is sensitive to the region with the relatively 
higher virtualities for which the soft photon can still resolve the soft and hard-collinear structures, i.e. 
where soft photon is sensitive to the rapidity gap arising due to the soft spectator scattering. 



6 Phenomenological analysis 



In order to make numerical analysis we exclude the SCET FF T\ using the expression in E q.(91 1. Consider 



for simplicity the amplitude F 3 . Using the soft spectator contributions given in Eq.(75l we obtain 

a v 



;F 3 (e,Q z )=g 3 (z,Q) 



-C 3 (s,t)J"i + &(xi) * H 3 (z,Q 2 , Xi,yi) * 



(132) 



where for the hard spectator scattering contribution we used expressions given in Eqs.(|72|). Substituting 

assume t 

T2(s',t) 



expression for T\ defined in Eq.(91| into Eqs.(132) (we assume the different value of s = s' for WACS) 
we obtain 



777/ 



53(2, Q) + - -C 3 {s,t)- 

7T S G 2 (s',t) 



+*(!/<) * <^ H 3 (z,Q 2 ; Xi,Vi) 



H 2 (s',t; Xj,yj) 



7T S 



-C 3 {s,t)1l{Q) + * H 3 {z, Q 2 ; x u Vi ) * 9( Xi ). 



(133) 



(134) 



In order to obtain the last equation we used the estimate ( 131 1 in which g 3 is zero, substituted 
T 2 (s' ,t)/C2{s',t) = TZ(s',t) ~ TZ(Q) and we also neglected the contributions of order 0{a 2 8 ) in the hard 
coefficient functions #2, 3, keeping in mind that the leading order contributions are H$ t M ~ 0{a s ) while 
H2 ~ 0{a 2 ) and C2 ~ 0(aP s ). An analogous expression can be also obtained for the TPE amplitude 



SG^ e and reads 



SGf/(e,Q 2 ) 



- G M {Q 2 ) In ^ 

IT ft* 



+ -C M (z,Q 2 ,f J L )TZ(Q) 



V(yi)*H M (z,Q 2 ; x u yi) * 



(135) 



where we used for g\ the expression from Eq.(131) and assume that Hm ~ 0(a s ). 

The formulas presented in Eqs.( 134|135 1 will be used for the numerical estimates presented below. 
We like to emphasize that these expressions are of course model dependent. We expect that the largest 
ambiguity arises from the amplitudes g± t3 which were computed in the hadronic model in Eqs. (1311, 
and which describe the TPE subprocess with one soft photon. 
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The other terms in Eqs.( 134|135 ) describe the contributions with the two hard photons and can be 
fixed more accurately. For the ratio 1Z we used numerical values given by the fit ( 107 1 with the parameters 
c = 0.10 ± 0.01 and A = 1.08 ± 0.03 GeV. In the numerical calculations we use for the scale ^ = 1 / 
that yields /i = 630 — 450 MeV for Q 2 = 2.5 — 5 GeV 2 . We observed that variations of the value of 
this scale by ±100 MeV do not provide any significant numerical effect. The explicit expressions for the 
coefficient functions Hm,3 and distribution amplitudes ^ can be found in Ref. 32 33 
be expressed in terms of one scalar DA for which we use the model 



The DA * can 



^(11,12,13) = /jv 120x1X2X3 {1 + r + (l - 3x 3 ) + r_(x 2 - Xi)} , 
where the low energy constants /jv, r+ and r_ have been estimated at /i = lGeV as [30] : 



f N = (5.0 ± 0.5) x 10~ 3 GeV 2 



0.35, r_ 



1.37. 



(136) 



(137) 



The scale of the running coupling a s (^) was set to be ~ 0.6 Q 2 and we use the two-loop running 
coupling with a a (p% = 1.5 GeV 2 ) = 0.360. 

In our numerical calculations we also use the ratio Ge/Gm which is fixed from the experimental 
data [l}{5] and in the numerical calculations we used the following fit of the data for Q 2 in the interval 
0.5 - 8.5 GeV 2 : 

Gl ^(l-ahV/AU), (138) 



R 



G 



M 



with a = 0.070 ± 0.007, A R = 0.54 ± 0.03 GeV. The other parameters which we need are 

m = 0.938 GeV, /i p = 2.7928, a = 1/137. 



(139) 



Finally let us inspect behavior of the amplitudes ( 134|135 ) in the limit s 1 (z(e — >• 1) — >• 0) where 
TPE corrections must be consistent with behavior Eq.(14|. In this limit using expressions for Cm, 3 given 
111 Eqs.( [66]67] we obtain 



6GT(e = l,< 



~F 3 (e^l,Q 2 



where we used that hard spectator contributions vanish separately at e = 1 

lim * H M ,3(z, Q 2 ; x h yi ) * *(x 4 ) = 0. 



(140) 



(141) 



This can be easily shown using the explicit expressions from Ref. 32 33 . Using Eqs.(140) one finds that 
Eq.(fT4j) are fulfilled. 



6.1 Reduced cross section 



The formula for the elastic reduced cross section presented in Eq. ( 20 ) can not be compared to the 



observable cross section because of IR-divergent amplitudes. For the real parts, the IR divergence arising 
from the elastic next-to-leading (NLO) QED diagrams is canceled by IR-divergent contributions arising 
from the interference of the bremsstrahlung and Born diagrams. Therefore one must include also the 
contribution of the inelastic diagrams. After that the reduced cross section can be presented as following 



exp I7 



(142) 



where the Born cross section reads 



17 _ 



= G 



M 



-G 



(143) 



The contribution of the elastic amplitudes in Eq.(142| is included in 5|L The analysis of the radiative 
corrections for the existing experimental data has been carried out using the Mo and Tsai (MT) formulas 
10 . Let us write in this case the reduced cross section as 



l7,MT 



rei.MT 
5 2 7 



rep,MT N 
brcms , 



(144) 
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In this case the extracted values of cr fl 7 ' 1V11 can be considered as experimental results. If one uses a 



different approach to compute the radiative corrections with a different elastic input then the extracted 
value of the Born cross section is different. This difference must be considered as effect related to the 
TPE correction in our case. In what follows we use the same expressions for the elastic contributions 
other than two-photon exchange, and for the bremsstrahlung contribution as in MT-formulas. Therefore 



from Eq.(142| and Eq.(144) we obtain 

4 7 ( 



I7.MT 

T R 



°2 7 °2 7 



ei.MT 



(l + 6% c - S% cMT 



(145) 



where the FFs Gm and Ge on the rhs must be considered as unknown quantities. Using Eq.(20l one 
can write the explicit expressions for the TPE contributions as 



a^S^ = 2G M BB 



scf; 



v ~ 
e^F 3 

Tfl 



2-G E Re 

T 



8GT 



v - 
— F 3 
to 



(146) 



It is convenient to present the analytical expression for the MT term 5 2 ^ C ' MT in the similar form. Using 
the analytical expression presented in [38] we can write 



4 7 



,MT 



2G M Re 



1 rtpe.MT 



2-G E Re 

T 



Ge k P 7 



!,MT 



The analytical expression for the factor 8 2 ^ / c ' MT reads 
Re(5 tpc,MT = 2a 





§ 




§ 




{'4'° 








- Li (D 




u 




u 





- i In 2 - + Re Li 
2 s 



iln 2 " 
2 u 



(147) 



(148) 



where A 2 is the soft photon mass (IR regulator), s = s — to 2 , u — u — to 2 and Li(z) is the Spence function 
(dilogarithm) defined by 

In(l-t) 



Li(z) 



dt- 



t 



(149) 



The result presented in Eq.(148) is obtained by computing the box diagrams with nucleon intermediate 

uswer 

(150) 



state in the soft limit 10 . However the expression in Eq.(148) is incomplete because the exact answer 

rtpo,MT 



has one more term 38 



4 p 7 e ' MT 



J 2 7 



How ever the factor air, which originates from the crossed box diagrarrj^J is absent in the formulas in 
Ref. [10] and therefore is excluded from the expression oiS^?' MT . 
Using Eqs.(146), (147) we can present the rhs in Eq.(145) as 



a R = <*R 



2G M Re 



sdf; 



r» X*pe,MT 



2 — G £ Re 

T 



6G^ C 



rn 



V Z7< n rtpe,MT 



(151) 



The amplitudes SG^, SG^ and 5^ C ' MT in Eq.(151) depend on the IR photon mass A 2 which however 
cancels in the difference, we discuss these subject m details in A ppendix D. 
The TPE contribution associated with the SG^ in Eq.(151) 



8GT 



^F 3 = SG\f - (1 + r)5F 2 + -F 3t 



(152) 



contains the unknown amplitude 8F2. This amplitude depend on the matrix elements of subleading SCET 
operators which we did not conside r in o ur analysis. We will assume that the corresponding difference 



sgT + ^f 3 -g e ^ m > 



therefore 



in Eq.(|15l| behaves similar to the FF Ge — Gm — (1 + t)F 2 = RGm and 

(153) 



5GT + — 2 F 3 - G4C' MT - -O(R). 

Tfl 2 ' 7T 



9 More explicitly, the factor aw in Eq. ( |l50[ l originates from analytic continuation of In 2 terms when evaluating crossed 
box diagram from the direct box diagram as shown in Ref. [49| 
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In this case the correction originating from such term in Eq.(151) is relatively small 



2-Ge Re 
r 



6gT + —^h - g e 6$; 



c,MT 



2-G M - x 0(R 2 ), 

T IT 



(154) 



because the ration R is small and we can neglect it. Therefore in our numerical analysis we use for the 
reduced cross section the following expression 



l7,MT _ n 2 
5 ( / 



•R 



M 



: G 2 E + 2G M Rc 



1 



SGf; + e^F 3 -G M ffi MT 



(155) 



The values c]j 7 ' MT must be considered as experimental data. Taking into account that the ratio 



R = Ge/Gm is also measured, see Eq.(138l, one finds that only the FF Gm is an unknown quantity on 



the rhs of the Eq.(155). Therefore using Eq.(155) we can extract the value Gm using the data from the 
kinematical region where our approximation is valid: s ~ —t ~ — u 3> A. We use the data from Refs. 7(9 . 
In Figs |13|14| we show the fit of the experimental data for different values of Q. The shaded area shows 
a transition region where the Mandelstam variable u is already quite small and one may expect sizable 
contributions originating from the higher order and, probably, power corrections. The corresponding 
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Figure 13: The fit of data for the ep — > ep elastic scattering reduced cross section from (7| for different 
values of Q 2 . For the fit we used only the data points which are lying on the right of the shaded area. 
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Figure 14: The fit of Gm using the data [9] with the computed TPE corrections. The black solid line 
shows the ub, with G'ff' extracted from the linear fit in 65 , black dotted curve represents the fit of the 



data points which are lying on the right of the shaded area. 



27 



values of | it | are shown at the bottom of the plot. For the fit we used only the points which are lying on 
the right of the shaded area. First we observe that the obtained TPE contribution is relatively small. It 
turns out that the soft spectator contribution by coincidence is numerically very close to the MT result 
for the considered range of Q 2 . 

As we can see from the plots the obtained behavior of <Jr(s) (solid black line) is linear in the region 
of e where the calculation can be applied and describes the difference in the slope between polarized and 
unpolarized data quite reasonably for SLAC data with large error bars and worse for the JLab data. 
However in the last case we have to take into account that, for instance, for Q 2 — 4.1GeV 2 the second 
point corresponds already to a relatively low values of u. 

The main trend is that in the region of small e the slope must change and this means that the 
reduced cross section can not be described by the linear function in e at large Q. Such nonlinearity is 



still relatively small at Q 2 = 2.64 GeV 2 (the linear fit is shown by the blue dashed line in Fig 14 ). But 
at Q 2 =4.1 GeV 2 this effect might already be sizable. Such a scenario implies a large corrections when 
u ; — > 1 GeV 2 , i.e. in the backward region. Theoretically it may be interesting to go beyond the formalism 
presented in this work and compute the TPE corrections using QCD factorization in the region with 
small \u\ < m 2 in order to have description for the whole region of e. 

The extracted values of Gm are shown in Table ljand for comparison we also provide the results of 
different phenomenological fits made in Refs. [65}|67 . 



Table 1: Results for magnetic form factor obtained in this work from the fit of SLAC data [7] (black 
solid line in Fig 13) and JLab data [9] (dotted curve in Fig 14). The extracted values denoted as G^ ct K 



For comparison we show the phenomenological extractions presented in Refs. 
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In Fig{l5]we show the TPE amplitudes which are defined as 

Y M = (sGf; - GmIs% cMT ) /G m , Y 3 = SF 3 /G M . (156) 



These ratios are plotted for three fixed values of Q 2 — 2.64, 3.2, 4.1 GeV 2 as a functions of e. For the FF 
Gm we take the values shown in Table [TJ As we can see these functions show weak dependence on Q 2 for 
given values of Q 2 . At Q 2 = 2.64 GeV 2 these functions has been estimated from the data using certain 



assumption about their behavior at e = 1, see e.g. 65 68,69 . Thus the present calculation shows that 



the assumption that the amplitudes Ym 3 vanish separately at e — > 1 is probably too strong. The absolute 
values of the Ym.3 in Fig |15| are also much smaller than the phenomenological ones shown in [65] . 

The other good observable to access the TPE corrections is to measure the ratio of the elastic e ± p 
cross sections 

d<jp+ n 

R± = (157) 

d(j e - p 

If the corresponding reduced cross section a 1 ^ extracted with MT radiative corrections (lOj are incomplete 
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Figure 15: TPE amplitudes (156) as a functions of e at fixed values of Q 



then this ratio will differ from unity due to the TPE contribution. In our case using Eq.(155) we find 
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(158) 

Using the values Gm obtained from the analysis of the JLab data we obtain the ratio R± shown in Fig[l6| 
as a function of e for three different values of Q 2 = 2.64,3.2 and 4.1GeV 2 , respectively. For comparison 



we also show the predictions made in 65 on the basis of linear fit for the same values of momentum 



transfer. Because we obtained a relative large difference between the slopes for the reduced cross section 



we also obtain a large difference in the estimate of the ratio (158). It is interesting that the absolute 



value of the obtained TPE correction within the factorization formalism changes very slowly with respect 
to Q 2 . 




Figure 16: The ratio R± for different values of Q 2 (in GeV 2 ): Q 2 = 2.64 (black solid curve), Q 2 = 3.2 
(blue dashed curve) and Q 2 = 4.1 ( red dotted curve). For comparison we also show the results for the 



linear fit obtained in 65 . The value of e are restricted to the range \u\ > 2GeV 2 . 



Let us stress that the above results obtained within the simple model of the SCET amplitudes gi 3, see 



Eq.(131). One cannot exclude that this rough model underestimates the contribution of the amplitudes 
gi^. Nevertheless the above calculations clearly show that the TPE contribution associated with two 
hard photons cannot describe the reduce cross section at large Q 2 being consistent with the assumption 
about the linear behavior of <7n(e) in the whole interval < e < 1 at large fixed Q 2 . Therefore, in 
order to understand the TPE contribution better it is very important to obtain a realistic estimate of the 
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amplitudes (71,3. It turns out that the important information about the amplitude gs can be obtained 
from the analysis of polarization observables. 



6.2 Recoil polarization observables 

The GEp-27 collaboration has measured the e-dependence of the recoil polarizations at a fixed value of 
Q 2 = 2.64 GeV 2 |6j. These data provide additional important constraints on the theoretical analysis of 
the TPE amplitudes. 

Let us to start the discussion with the longitudinal polarization Pi which in the elastic approximation is 



given in Eq.(22). Taking in to account the inelastic contributions and assuming that radiative corrections 



for the experimental data have been performed according to MT formulas we can follow the same line 
as for the reduced cross section. Introducing the new FFs Gm and Ge and expanding the reduced cross 
section we obtain 
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(159) 



Using Eqs.( l~4|l40 1 one can easily find that in the forward limit e — > 1 this asymmetry does not vanish 
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(160) 



Taking into account the relati ve sm allness of 0(R )-contributions and our assumption (154) we can 



conclude that the ratio in Eq.(159) is dominated by the first term 
approximated as 

pi j, Bom ■ - - 7 '■ 
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F3 in the brackets and can be 



(161) 



Taking into account that F3 is IR-finite and F^(Q 2 , e) — in the MT-calculation we can u se P in order 
to extract information about the amplitude F3. The kinematical factor e 2 appearing in Eq.(161) reduces 



the possibility of s uch analysis only to the region of relatively large values of e where our approximation 



is valid. In Fig 17 we demonstrate our estimates for the Pi/P®° rn as a function of e for different values 
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Figure 17: The ratio Pi/P l Born as a function of e for different values of Q (left) and comparison of 
our estimates with experimental data [6| at Q 2 = 2.64GeV 2 . The shaded area around the dashed line 
demonstrates the ambiguity arising from T\ . The solid line shows the ratio where T\ obtained from the 



fit with the kinematical power corrections, as in Fig 11 



of Q 2 . It turns out that the hard spectator contribution F^ is quite small (red dotted line in left plot 
) compared to the soft spectator contribution F^. From Figjl7]we also see that the dependence on Q 
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is rather weak. However this function has a steep behavior in e in the vicinity of e = 1. Such behavior 
originates from the e-dcpendence of the hard coefficient function 63(2), see Eq.(67). On right plot in 
Fig 17 we show a comparison of our estimate at Q 2 = 2.64GeV 2 with the existing experimental data [6j. 
For comparison we show the solid and dashed curves corresponding to the two different fits for J-±: with 
and without the kinematical power corrections, respespectively. The shaded area around the dashed 
curve shows la error bands arising from the ambiguity in extraction of T\ from the WACS data. As 
we discussed in Sec|4]the kinematical corrections increase the extracted value J-\ and therefore the solid 
curve lies a bit higher. The calculation shows the qualitative trend of the data, which increase when 
e — > 1. However our estimates are considerably smaller then the data point at large e. This observation 
can be considered as an indication that the contribution associated with the amplitude 33 is not small 
and probably even quite large. If such scenario is realized then the similar contribution to g\ must also 
be large in order to describe correctly the reduced cross section. 

The other important issue is that the TPE amplitude ^2 ^3 does not vanish at e = 1 as it has been 
assumed in phenomenological analysis in Ref. [65]. Notice that this behavior is closely associated with 
soft spectator contribution. 

Taking into account the effect of the computed TPE corrections in the transverse polarization yields 
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In the limit e — > 1, the TPE corrections vanish: 
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(162) 
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Taking into account assumption (154 1 and the fact that the second term in Eq.(163l is proportional to 



R we obtain the following estimate for the TPE correction 
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(165) 



This allows one to conclude that the expected correction can be smaller than 1%. The results of GEp- 
27 collaboration [6j are in agreement with this conclusion. Within their error bars of order 1% this 
experiment does not see any systematic TPE effect on this observable. 



For illu stration in Fig 18 we show the effect provided by the known TPE contribution on the rhs 
of Eq.{l63) jf-^h + SG M ~ G M ^7 C ' MT - The solid line is an ^independent fit of R. The dotted 
line shows the effect of the computed part of the TPE correction, the curve is restricted by region 
e > 0.6 (\u\ > 1.83GeV 2 ) We can clearly see that this effect is smaller than the experimental error bars. 



7 Conclusions 

In conclusion, we have discussed the calculation of the TPE amplitudes using the QCD factorization 
approach in the region where all Mandelstam variables are large (wide-angle scattering) s ~ — t ~ 
— u ^> A 2 . The leading power behavior of these amplitudes are described by the sum of the hard 
and soft spectator scattering contributions. The hard spectator contribution at leading order accuracy 



in a s was studied before in 32 33 . The soft spectator contribution was estimated before using the 
handbag approach in 38 . The aim of this work is to develop a more systematic description of the latter 
contribution in order to reduce the model dependent assumptions. 

In our work we used the SCET framework and suggested a more rigorous factorization formula for the 
description of the soft spectator contribution. We carried out the leading power analysis in the hard scale 
1/Q 2 for the TPE amplitudes and obtained that the soft spectator contributions to the TPE amplitudes 
can be presented as a sum of two terms associated with the different hard configurations of the photons. 

If both photons are hard then the corresponding amplitude is described by the hard coefficient function 
and SCET FF denoted by J-\. This function describes the soft overlap contribution between incoming and 
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Figure 18: The ratio — /z p ^/t(1 + e)/2sP t /Pi as function of e for Q 2 = 2.5 GcV 2 . The data are from 
GEp-I (blue triangle) [l][2], and GEp-27 (red circles) [6]. The solid line is an e- independent fit, the dotted 
(green) line shows the effect of the known part of the TPE-correction. 



outgoing proton and depends on the hard-collincar scale of order AQ. Therefore this function depends 
only on the momentum transfer Q. The £-dependence is completely described by the hard coefficient 
function and can be computed in perturbation theory. 

If one of the photons is soft then the corresponding SCET amplitudes denoted by gi^ are different. 
These functions describe complicated low-energy hard-collinear and soft dynamics. In case of the soft 
spectator scattering the soft photon is involved in interactions with the hard-collinear and soft constituents 
of the proton. In this case the e dependence arises due to these long distance interactions. 

At moderate values of Q 2 when the hard-collinear scale is not large QA < m 2 the SCET amplitudes 
describe the long distance interaction and cannot be computed in pQCD. In order to estimate these 
quantities we used different approaches. 

The contribution of the FF T\ can be estimated using the universality of its definition in the SCET. 
It turns out that this FF also arises in the factorization of the wide-angle Compton scattering amplitude 
(WACS). We considered this process at leading order in 1/s and in the coupling a s and provided the 
SCET factorization for the soft spectator contribution of the WACS amplitudes. We found that the 
soft spectator contribution of all three dominant amplitudes are described by the same SCET FF T\. 
This allows us to establish the specific relations between these independent amplitudes. Assuming the 
dominance of the soft-overlap mechanism we extracted the values of the dominant amplitudes in the 
region Q 2 = 2.5 — 6.5 GeV 2 using the JLab Hall-A data for the WACS cross section. The contribution 
from the SCET FF T\ can be related to the physical WACS amplitude and this provides an estimate of 
this part of the TPE contribution. 

A description of the amplitudes g 13 is more challenging because the corresponding configuration is 
very specific and it is difficult to find a process where they also enter. Moreover, the amplitude gi also 
includes the QED IR-divergent part which must cancel against inelastic terms. Therefore in order to 
estimate these amplitudes we assumed that the dominant contribution might be associated with the 
exchange of the very soft photon momentum. Such photon interacts with the total charge of the proton 
and cannot resolve the charges of the constituent quarks. In this case one can estimate the 171,3 using an 
effective theory with nucleon degrees of freedom similar to the hadronic model which has been used for 



TPE calculations at relatively small Q 12 17 . We suggested a somewhat similar calculation using for 



simplicity only the nucleon as intermediate state. The resulting expression for 171 provides the correct 
QED IR-behavior and does not include any large logarithms which can be associated with the contribution 
of the hard region (in contrast to the MT calculation flO]). This point is very important because the 
subprocesses with the hard virtualities are described by the hard coefficient functions. 

The amplitude 53 does not have any IR-divergence and we obtained that g^ ~ in the framework of 
the simple hadronic model with the nucleon as intermediate state. This allows one to conclude that the 
dominant contribution to the amplitude 53 may arise from the region where the soft photon virtuality 
is of order A 2 . In this case the soft photon can still resolve the hard-collinear and soft structure of the 
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proton arising from the soft spectator scattering. An investigation of this dynamics is definitely beyond 
the simple hadronic model used here. In our work we do not investigate this domain and accept as a first 
approximation the results of the elastic hadronic calculation. 

Summarizing, we defined the two dominant TPE amplitudes in the framework of the QCD factoriza- 
tion approach. We found that the contribution associated with the two hard photons can be computed 
in terms of reasonably known nonperturbative quantities: nucleon distribution amplitudes (hard specta- 
tor contribution) and SCET FF T\ (soft spectator contribution ) which can be fixed from the WACS. 
The largest ambiguity arises from the asymmetrical configuration with hard and soft photons which is 
described in terms of the SCET amplitudes 31.3. The latter can only be computed in the framework of 
low-energy models. For simplicity, we consider a scenario when these amplitudes are dominated by a 
very soft photon which already interacts with the hadron as a point-like particle. The computed TPE 
corrections have been used in a phcnomcnological analysis of the experimental data in the region where 

1 1*| > 2.5 GeV 2 . We obtained that the corresponding TPE correction are linear in e but relatively 
small in magnitude so that already for Q 2 > 2.5 — 3 GeV 2 a description of the reduced cross section is 
not compatible anymore with the pure linear fit used for the Rosenbluth separation in the whole interval 
< e < 1. Qualitatively this observation coincides with the hadronic model calculations which also show 
the similar trend for the moderate values of Q 2 . 

We investigated also the ratio R e + / e - which is sensitive to the same combination of the TPE ampli- 
tudes as in the reduced cross section. We obtained that our predictions with 53 ~ for the R e + i e - at 
Q 2 = 2.5 — 4.1 GeV 2 are considerably smaller than the estimates originating from the linear fit of the 



reduced cross section, see e.g. 65 



However in the considered model such a small TPE effect in or can also be a consequence of the 
assumptions about the nonperturbative amplitudes g\^- We suppose that the used hadronic model can 
underestimate the contribution of these functions. The important information about the amplitude 53 
can be obtained from the polarization observable P\. We found that the ratio P;/P ; Born is dominated 
by the soft spectator contribution. However using (73 ~ our estimates yield a value for the ratio of 
0.5 - 1% for Q 2 = 2.5 GeV 2 and e > 0.7, falling short of the data. Further measurements of Pi/P l Born 
for different values of e and Q 2 with better accuracy may help us considerably to constrain the unknown 
173 and reduce the uncertainty in our analysis. 

The small effect of the TPE corrections in the transverse polarization P t can be explained by the 
suppression of the helicity-flip amplitudes comparing to helicity conserving ones at moderate and large 
Q 2 values. This assumption is realized for the proton FFs in the smallness of the ratio Ge/Gm and we 
expect that such suppression also holds for the corresponding combinations of the TPE amplitudes. 

Our work shows that even for the wide-angle kinematics where s ~ — t ~ — u ^> A 2 the hard TPE 
contribution ( both photons are hard) can not provide the dominant effect. We expect that the relatively 
large contribution also arises from the region where one of the photons is soft and associated with the 
nonperturbative hard-collinear and soft dynamics. Such scenario is supported by the large TPE effect 
in the longitudinal polarization Pi which is considerably larger than predictions obtained using only the 
hard TPE contribution. The other indication is the relatively small hard TPE contribution to the reduced 
cross section in the large e region. Further measurements of the e dependence of the reduced cross section 
or the ratio R± at intermediate Q 2 values will test any non-linearities due to TPE effect. 
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8 Appendix A 



In order to obtain the expressions ( 39 ) for the FFs it is enough to consider the matrix element of the 
electromagnetic current. We also use 

P = Q^ + ~+0(m 2 /Q 2 h P ' = ^ + Q^ + 0(m 2 IQ 2 ). (166) 



33 



Lets us define 

N(p>) = N(p')^ + N{p')^ = N' + + N'_, 
N(p) = ^N(p) + i^N{p) ^N + + N_. 
Then from the equations of motion for nucleon spinors one obtains 



that gives 



Therefore 
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N(p')lN(p) = N' + 1N + + 0{m 2 /Q 2 ). 
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On the other hand the leading ordder SCET expression gives 
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Comparing (172) and (1731 we obtain 
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(167) 
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(169) 
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where F 2 denotes the contribution from the subleading SCET operators. 



9 Appendix B 

The combination of the matrix elements which define the SCET amplitude g\ reads 

(p'\ T{O^t(0)S o (0)} b) SCET + (p'\ T{O»Yj(0)Sz(0)} \p) SCET = a g q 1 (z, Q,fi F ), (175) 



where the second term (p' \ . . . \p) denotes the antiquark contribution. We consider only calculation of the 
fist matrix element on lhs ( 175[ ). The explicit expressions for the diagrams in Fig[5]can be written as 
follows 



ie 2 e 2 u 2e - 



(27T) 



2„2 r ,2e 



D*„. + D 



2 b 




(2ir) D s " ' s " J [I 2 -X 2 + ie] [(In) + ie] { [(vl) + ie] [- (vl) + ie] 
Here we assume that the renormalization scale p 2 is defined in MS-scheme 



/2 2 = 4e-^«- ln47r . 



(176) 
(177) 

(178) 



One can easily see that each integral Dik represents the evaluation of the T-product of the WLs associated 
with the couple of the light-cone vectors: v and n, v and n and so on. In fact, soft photons can be 
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decoupled from the hard-collinear quarks performing the similar redefinition of the quark fields as we did 
for the lepton fields in Eq.(49). Then one obtains that the diagrams in Fig[5] represents the calculation 
of the T-product of the light-like WLs. 

Computation of the integrals in Eqs. (176) and (1771 are very similar therefore it is enough to 
consider only one combination. Each integral Dij are not well defined even in dimensional regularization 
and therefore one needs to introduce some additional regularization. Consider the integral in (1761 and 
let us introduce the following regularized expressions 



D 



ith 



Jla 



Jib 



la,b 



dPl 



ie 2 e 2 u 2e 



£n7 M &i Jla,b 



1 



(vn) 



A 2 + ie] [(In) 
d D l 



t_ + ie\ 
1 



(v I) - r + + ie] 
(vn) 



(179) 



(180) 



(181) 



[P — A 2 + ie] [(In) - r_ + ie] [(vl) - r+ + ie] ' 

where t± are additional regulators required in the intermediate calculations. With these regulators we 
can compute the integrals separately. Consider Jib, using light-cone vectors fi and v as basic light-cone 
vectors: 



1= (ln)-— + (lv)- 
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we rewrite the integrals in terms of these Sudakov variables 

1 1 



(vn) 



J 



where we used that 



16 



d°l 



dl+dl-dl i 



[zl+l- A 2 ] [l_ - T-] [1+ - t+] ' 



1 



-dl+dl-d 



D-2 



l±, (vn) 



2(k'p) 



-u 
y2 



(182) 



(183) 



(184) 



(vn)-^-"— v -' v Q 2 Q 2 

The obtained integral can be easely computed using residues and then integrating over transverse 



momenta: 
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The second integral can be computed in the same way, we need only to take into account that 
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where we keep only the real part. Then 
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The second sum can be computed similarly and the sum of the all terms reads 

D la + D lb + D 2a + A26 = ~ Ll^n e\ | \ In z + In z In ^| } ( 1 !tt ) 

The antiquark can be computed in the same way. The UV-pole 1/e is removed by the renormalization 
and the final expression for the SCET ampltude g\ reads 

gl(z,Q,[i F )=\nz\n^F?. (194) 

10 Appendix C 

In order to see the cancellation of contributions with the hard-collinear photons we suggest to compute 
the diagrams in Fig|4] using the strategy of regions 70 71 . Using the IR-regularization with off-shell 
external momenta we obtain 

Di,2 = -i-^jjn J ^ ( 195 ) 

where subscript 1,2 denotes the box and the crossed box diagram in Fig|4j respectively. Corresponding 
integrals read 

f UvTil^yfu, i (/ + jj) 7%-> nqm 

1 J [p-2(ik)-ki] [i 2 m+p-p') 2 \ [i 2 + 2(i P )-pi\ ' 1 ' 

J = frll + f)rUv 1 ^ (/ + $ nq 7 ) 

' 2 J [P + 2 (ik')-ki) [p][{p> -p-iy\ [p + 2(i P )-pi] ■ [ ' 

where we assume that poles in the brackets [...] are always defined with +ie prescription. 
In the collinear to p region we obtain 



Jlcp - U^U V £ nl "tn / dl r _ 2{lk) _ ^ j { _ 2pl{l + p)] | [p] [p | 2( ^ } _ p 2 ] f" > < 1!,S 1 



and 

/" 2fc ,Al 1 f 2(Z+p) /i 

J 2 op=U^U v ^^nJ dl [2{w) _ ^ h2p , (?+p)] [Z2 +2 (/p)_^_] |- ,l!,<n 



Therefore 



1 f 2(/+p)^ 



[-2p'(Z+p)] l[J a ][l 2 + 2(Jp)-pi] 
2fc^ 2A^ 



[-2(^fc)-fc2] [2(;fc')-fci; 



(200) 



Consider expression in (...). 



k^n^ k_n^ 
( " ~ [-/+*:_ - k\] + [l + k'_ - kY] = H+ - k 2 Jk_] + [l + - k 2 Jk'_] ' 

Notice that the IR-regulator k\ introduces the small difference. Consider the soft limit 1^ ~ A in the 



collinear integrals in Eqs.(200). We obtain 



1 f 2p M 1 

( Jlcp + J2cp) s = Uyj v u v CnT^fi / * ^ | jjaj -pjj ) ( 202 ) 

h + - Ai/*-] + [i + wm J ■ (203) 
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This expression coincides with the soft limit taken in the original integrals Ji.2- Therefore in order to 

yields 

(204) 



obtain the collinear contribution we must subtract the overlapping soft part, see e.g. 59 60 . This yields 

\J\cp J2cp) {J\cp '^2cp) s — Hy'f u v £nT Cn Jcp: 



with 



Jcp 



dl 



1 



1 



2pf 



[-2p'(l + p)] [P] [P + 2(lp) - p\] [P] [2{lp) - p\] 



k 2 Jk'_) 



The integral J cp is not singular in the region of small l + and we can put k\ = that yileds 



[-1+] 



= -2m d{l+)n fi . 



Hence 



(205) 
(206) 

(207) 

(208) 
(209) 

[-P ± + 1-p + -pI}[1-p + -pIY (210) 

This integral is trivial because the poles in Z_ lie on the same complex semi-plane. Thus we demonstrated 
that after soft subtraction the collinear to p integral is zero. The same manipulations also can be done 
for the other collinear regions. 



Jcp — (— 4) 7T? / dl 



5(1, 



P 2 ] V 1 

dl_dl± 



1 (1++P+) 1 P+ 

-p'_(l++ P+ )] [P + 2{lp)~p\] f[2(lp)-pl) 



1 



1 



1 



Hi] V HI + J-p-4 



■Pi] 



[i-p-i 



■p\] 



= ~^P+ l dl-dl± 



1 



11 Appendix D 



The dependence from the IR-regulator A in amplitudes £G^° £ arises only in the SCET amplitude gi 



Using Eqs.(20), (74) and (131) we obtain 



SG^ e = ~gi(z,fio) + ... = -G M ^2 ln 

7T 7T jJ, Q 



(211) 



where dots denote the A-independent terms. Therefore for the corresponding difference we obtain 

2 "2 7 



Re[6G\ p ; - G M \8% cMT ] = -G M In ^ ln 

" 7T // 



Q; A 
— Gm In — ln 

7T S 



— Gm In Q - ln 
7T ^ 



(212) 



The expressions ( 39 ) and ( 76 ) for the FF F-2 and amplitude <5F 2 are incomplete because we neglected 



the contributions with the SCET subleading operators; Therefore we can consider cancellation of the 

IR-regulator A in the combination SGe — Ge^i^ in Eq.(155) only between the corresponding kinematical 



power corrections. Using Eqs.(76) yields 

a 1 



5Gf = SG M - (1 + t)SF^ ~ --- { gi (z,fio) + C m {z,im>)Fi + -C 3 (z)Fi\ - -C 3 (z)Fi, (213) 

7T r i. s > s 



Using Eqs.(20) and (131) we obtain 

a 1 



SG 



-- { 9i{z,Ho) + C M (z, ho)^i + -C 3 (2;)J"i) 

7T T L S J 



al , . a 1 „ , „, . A 2 
5i(z,^ )+ ••• = GA/((y)ln-o ln 

7T T 7T T /ig 



(214) 
(215) 
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Then taking into account (39) where we also neglected the subleading SCET operators we obtain 

1 „ 1 



Ge — G 



M 



Taking into account Eq.(148) yileds 



^ 1 rtpe.MTi 

Ge 2 d 2 7 J 



(1 + t)F 2 ~G m -(1+t)-G m 

T 



Gm In — ^ In 

7T r ^ 



-G 



Af ■ 



s 




1 \ A 2 


S 




Gm ) In — In 




u 




T J S 


u 



Gm In — «• In 

7TT /ig 



(216) 

(217) 
(218) 



From this calculation we can clarify why we keep unexpanded the argument of the logarithm in Eq.( 131 V. 
we need this otherwise the cancellation of the IR-mass A 2 will be incomplete because the expression for 
$2 7 in Eq.(148) includes such power suppressed contributions. 
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